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I. INTRODUCTION 


This final report documents four computer programs 
developed and delivered under NASA ARC contract #NAS2-9754. All 
of these programs provide copmputational assistance in the 
analysis of man/machine systems. The programs are as follows: 

(1) Modified Transfer Function Program (TF) 

(2) Time Varying System Response Program. (TVSR) 

(3) Optimal Simulation Program (TVOPT) 

(4) Linear Identification Program (SCIDNT) 

All of these programs have been installed and are opera- 
tional on the TSS-360 System at NASA-ARC. The Modified Transfer 
Function Program is a result of modifications to a previously 
delivered computer program described in ’’Human Operator 
Controller Model — (HOCOM) ; User’s Guide," by J. Mohr and 
K. Kessler dated April 1978. The program converts the time 
domain state variable system representative to frequency domain 
transfer function system representation. This program is 
described in more detail in Section II. 

The Time Varying System Response Program (TVSR) computes 
time histories of the input/output responses of the human 
operator model. It does this by means of a closed loop simula- 
tion derived from the system matrices calculated from the 
modified transfer function program. The link between these two 
programs is provided by writing the output of TF to output di$c| 
for storage. The simulation program reads this data and 
computes time histories. The simulation program was coded in 
such a way as to receive generic input data. That is, any 
linear system description can be written to disc (from card 
input, if desired) and read by the simulation program for 
subsequent processing. Specifically, Appendix A lists a stand- 


alone program which can be used to create a general output disc 
file to be read as input for the TVSR program. Due to the fact 
that the input data is obtained from a disc file, the TVSR 
program features a powerful mechanism to accurately simulate 
time varying manual systems. System matrices are read from disc 
for start and stop times. At appropriate intervals, the TVSR 
program interpolates between times to produce a smooth transition 
that approximates a time varying system. Several of these 
piece-wise linear constant coefficient simulations can be 
"pieced" together to create extremely versatile code. Section II 
describes the use of this linear simulation program (TVSR) . 

The third com.puter program is an optimal simulation program 
(TVOPT) , described in Section IV. It is similar to the program 
TVSR in that it produces time histories of system states 
associated with an operator in the loop system. The first major 
difference of these programs is that TVOPT explicitly accounts 
for the reaction time of the operator whereas the TVSR program 
approximates the delay/prediction process by a lumped linear 
system. The second basic difference is that TVOPT v/as designed 
to operate as a stand-alone program. The system description is 
directly read from input cards whereas the TVSR program obtains 
its data from disc data created by the transfer function program 
or by some other means (see Appendix A) . The third difference 
is that TVOPT cannot interpolate and, thus, represent systems 
constrained to a constant coefficient system. The last 
difference between these two codes is TVSR computes the expected 
value (and la uncertainties) for the states as well as sample 
time histories, whereas TVOPT is limited to only sample time 
histories. 

Both programs TVSR and TVOPT write the input/output time 
histories for the operator to disc for subsequent processing by 
the- fourth computer program, the linear identification algorithm 
(SCIDNT) . This last program is an open loop identification code 




which operates on the simulated data from TVOPT (or TVSR) or 
real operator data stored from motion simulators in an appro- 
priate fashion. This code produces estim.ates of the various 
parameters associated with the operator model and is a modified 
version of the SCI Identification Program (also named SCIDNT) 
that has been delivered to various other government agencies. 
This program is described in Section V. 

I Each program description is self-contained within the 
appropriate section. The description includes a sample input, 
the correct procedure to execute the various programs on the TSS 
system, library routines and input/output disc files where 
appropriate. In the descriptions to follow, it is assumed the 
reader is familiar Avith the theoretical considerations involved 
with the algorithms; in particular, it is assumed the reader has 
knowledge of the optimal control human operator model, which 
these programs utilize. 

In most cases, the notation used in the computer programs 
(as well as this user's manual) conforms to the widely used F, 

G, H etc. matrices defined in the optimal control literature. 

The operator notation also generally follows that used in the 
literature. The one exception is the SCIDNT program which uses 
different notation internal to the program. However, even in 
this case, every attempt \^a.s made to modify the output labels in 
this program to conform to those used in the transfer function 
and simulation programs. 


n 
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II. MODIFIED TRANSFER FUNCTION PROGRAM (TF) 


2.1 DESCRIPTION 

The Modified Transfer Function Program (TF) converts a time 
domain state variable system representation to frequency domain 
transfer function system representation. The transfer equations 
may be computed for: 

(1) plant/noise model only (no human model) 

(2) human model only - open loop 

(3) combined plant/noise/human model - open loop 

(4) combined closed loop plant/noise/human 

In addition, various frequency domain analysis aids are 
also computed. These include; 

(1) plant/human operator state dynamics eigensystem 

(2,) printer plots/tables of Bode magnitude and phase 

diagrams of any specified transfer function element 
arising from systems 1,2, 3, 4, above 

(3) poles/ zeros of any specified transfer function element 

(4) residue at all poles for any specified transfer 
function element. 

The TF computer code also computes the covariance matrix 
associated with the closed loop system dynamics (plant/noise/ 
human model).. 

The computer code has also been designed in conjunction 
with additional computer code (also delivered to Ames Research 
Center) called Time Varying System Response Program (TVSR) 
described in the next section. An option is available to write 
to disc, intermediate output from TF. The program TVSR then 

•Page blank 

5 - .C.:, 


PRECBBING 


reads this data and computes the time varying system response of 
user requested states and their associated la uncertainties. 

2.1.1 Program Functions 

The basic options are determined by the system configuration 
specified by the user. Transfer functions for multi-input/multi- 
output plant/noise/human models may be computed for the following 
cases (see Figure 2.1): 

(1) open loop plant dynamics only (open loop) 

(2) human operator dynamics only; (open loop) 

(3) human operator/plant dynamics (open loop) 

(4) closed loop system dynamics 

(5) covariance matrix 

The equations representing the above user options are given 
in Table 2 . 1 . • 

The basic structure of the program allows for repeated 
application of generic transfer code which utilizes a constant 
coefficient differential system of the form: 

X = Fx + Gu 

y = Hx + Du (2.1) 

where x is an n component state vector, u is an m compo-j 

nent input vector, and y is a p component output vector. 

The frequency domain representation of the same system in 
terms of the complex frequency s is : 

yCs) = T(s) u(s) (2.2) 
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where T(s) is a matrix of complex, rational expression, 
transfer functions between elements of the input u(s) and the 
output y(s). T(s) can be expressed in terms of F, G, H, and 
D as : 

T(s) = H(sl - F)’^G + D 


rriT 


[H adj Csl - F)G + A(s)Dl 


A(s) is the characteristic polynomial of the system and can be 
expressed as: 


A(s) = a^ + a^s + a 2 S^ + ... + 


or as 


ACs) = (s-p^) CS-P 2 ) . . . Cs-Pj^) 

The Pj^ k=l,...,n are the poles of the system, Eq. (2.1) 
The adjoint of (sI-F) is a matrix such that: 


ad j(sI“F) _ (•~T p")-! 
" aTsI 


In general, the elements of adj(sI-F) are polynomials in s 
and can be expressed in the form: 


[adjCsI-F)lij = b..^ + b.j^s ^ b..2^“ * b..jjS- 




A (s) 
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The Z--V the zeros of the transfer function between y. 

X J K 1 

and Uj . The number of zeros, !L, is less than or equal to n.. 
Given the matrices F, G, H, and D, the program can compute 


the n coefficients 


in the characteristic equation and the 


n+1 pxm arrays of coefficients b 


ijk' 


This generic transfer function code utilizes the Leverier 
algorithm described in the literature. Typically, systems on 
the order of 18-20 states can be determined with excellent 
accuracy utilizing both forward and backward computations. 
Numerical problems can and do occur \vith systems larger than 20 
states . 


In principal, then, transfer functions associated v/ith (1) 
plant/noise model, (2) human operator model, and (3) closed loop 
system are restructured so as to appear as Eq. (2.1). A slight 
variation on this common procedure is used to compute the open 
loop plant/noise/human model transfer functions. For this case, 
the open loop plant/noise transfer functions are computed 
followed by the open loop human model transfer function. The 
associated adjoint matrices are multiplied together to produce 
the resultant system (the characteristic equation for each is 
similarly multiplied) . 


2.1.2 Computational Aids 

Other computational aids for systems analysis are also 
available. These fall in the general categories of: 

(1) State Dynamics Eigensystem 

(2) Frequency Domain Analysis Aids 

(3) Covariance Matrix Calculation of the Closed Loop System 

(4) Optimal Cost Function of the Closed Loop System 
These are briefly described in the following subsections. 
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2. 1.2.1 State Dynamics Eigensystem 


Given the state dynamics matrix of the time domain repre- 
sentation of a linear system, the program is capable of 
determining the complete eigensystem. The method of Householder 
orthogonal transformations is used. The following parameters 
are computed: 

Cl) Eigenvalue (pole) locations in the complex plane 

(2) Natural frequency of each pole (magnitude of the 
eigenvalue) 

(3) Damping factor of each pole -RE(A)/|1 A |j 

(4) Eigenvector corresponding to each eigenvalue. The 
eigenvector is normalized so that its largest component 
has unity magnitude 

(5) Magnitude of each (possibly complex) component of each 
eigenvector 

(6) Phase angle of each component of each eigenvector 
relative to the largest component of its own eigen- 
vector 

2. 1.2. 2 Frequency Domain Analysis Aids 

The user may specify elements bf the transfer function 
matrix which he wishes to be analyzed in detail. For each 
element specified, the following options are available. 

Bode Plots 

The program can produce printer plots of the Bode magnitude 
and phase diagrams for the transfer function element. The user 
must specify the maximum and minimum frequency for the plot. 

The user also has the option to plot the transfer function 
element multiplied by an arbitrary numerator and denominator 
polynomial in s. 
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Zeros 


The program can compute the zeros of the transfer function 
element. These are the possibly complex roots of (s) = 0. 

The D.C. gains are also computed for convenience. 

Residues 

The program can compute the residues of the trahsfei* func- 
tion element evaluated at all of the system poles. The residue 
of the i,j element at the k th pole is defined as: 

^ijk ' 

s * Pk 

2. 1.2. 3 Covariance Matrix Calculation of the Closed Loop System 

The program computes the covariance matrix associated with 
the closed loop system. It utilizes the Q-R algorithm 
developed as part of the OPTSYS computer program. This algorithm 
fails with, zero eigenvalues. A standard procedure is to introduce 
small perturbations on the order of 10"^° to 10"^® (depending on 
the size of the matrix) in selected elements of the system matrix 
so as to slightly move the eigenvalues away from zero. 

2 . 1 . 2 . 4. , Optimal Cost Function of the Closed Loop System 

The program computes the individual costs JU, JX as well 
as the sum J (as defined in Table 2.1). 

2.2 UTILIZATION 

The Modified Transfer Function Program (TF) is very similar 
to a previous delivered transfer function computer program (see 
Section I) . Basically, this new version replaces the augmented 
states, labeled x^ in Table 2.1 of the above referenced document. 
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by new augmented states arising from a linear system approxima- 
tion (Fade approximation) for the transfer function between the 
control output, u, and the desired operator output, u^. To 
see how this is accomplished, refer to Figure 2.2. Here, the 
Standard Optimal Control Model (SOCM) is shown together with the 
linear system approximation required for the TF program. These 
blocks shovm in Figure 2.2 must be lumped into an equivalent 
system; given by the augmented states listed in Table 2.1. 

Combining the filter states with the augmented states 
yields the open loop transfer function of the operator in the 
form: 


where 


d 

3t 


A 


X 


X 


a 


u 

= 


f \ 


F (S-rah)H^ 


X 


K 

0 


y 


C-DaA . HJ 


X 

X 


+ BC V, 


F = Fj^ - KH^ - CGj,-KD^)D,A 


The open loop operator transfer function, u/y, is now 
computed rather than that discussed in the above referenced 
report. The human operator/plant dynamics (open loop) option 
now relies on u/y (given above) in the calculations. 


Example 

An example is given to show how the input matrices (provided 
by the user) are determined. 
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REQUIRED FOR TF • 

MOTOR 

V NOISE; V, 



APPROXIMATED BY LINEAR. 
LUMPED. ARBITRARY ORDER, 
TRANSFER FUNCTION WITH 
STATES X, 

a 


Figure 2. 2 Human Operator Models 
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plant 


Xi = -2X^ + w 
Xj = Xj * u 

y = X2 * Vy 

Yp (perceived) 


; =0 


’ V 


■= 0 


= y(t-T) ; T = delay time 


The user must input values for 


W = E (ww^ ) 


Vy = E(Vy Vy‘) 


as well as Q, g arising from the optimal cost 


J = j^llXllQ + llul|g}dt 


To obtain a linear system approximation to the predictor 
(in the Kalman filter) and the delay in the control, we approxi 

mate the predictor as 1+Ts and the delay from a second order 

Fade approximation as 


f(T) s 


6 - 2Ts 

6 + 4Ts + T^s^ 


The predictor and Fade approximation are lumped into an equiva- 
lent transfer function of the form: 



6 - 2Ts 

_6 + 4Ts + T^s^ 


(1 + Ts) 


6 4Ts - 2T^s^ 

6 + 4Ts + T^s^ ORIGINAL PAGE IS 

OF POOR QUALTIY 
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and BC is selected to produce "motor noise," v^, of suffi- 
i cient amount desired by the user. 

i ' Observe that the prediction and delay functions of the 

I : ' operator are lumped into one transfer function, currently up to 

I order 5. This is only possible because of the linear assumptions 

i 

I for these processes. 

■ ■ I The output data from the modified transfer function program 

j is identical to that used in the original version. The input 

i 

i data is slightly modified to include data from the augmented 
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states given by the matrices , G^, r^, D^, and 

(= cov ■'^A^ • 

Features 

There are several features embodied in the code that require 
some discussion. 

(1) It is possible to iterate on the values of and 

so as to converge to a solution (see Table 2.2, Group 1) 

Vy = Py E(y^) 

and 

\ ' Pa 

Either V (only) or V (only) or both, may be iterated. The 

y ^ ^ ^ 

iteration utilizes the subroutine MITER which in turn calls 
subroutines GRADE, ITERCK and AITERY . The algorithm for the 
iteration (in GRADE) for both and is 

V(new) = 10**[logj^Q V + log^o V(old)] 

where V is a generic representation of each diagonal element of 

✓V 

the covariance matrices V or V . V is given by 

/ * 


XCOV (i,i) ; i = 1,2,. ..,n 
n = # of elements in 
= constant 

XCOV = ECy^ yj); (see Table 2.1) 
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a 

V. = XCOV Ci.i) ; i = 

= constant 

ci 

j = # of measurements of operator (size of Hj^) 


m-j = # of elements in 

a 


Modifications to this algorithm may be included in GRADE to 
improve convergence. 

The "stopping" criterion (determined in ITERCK) is to check 
the diagonal elements for 


V(new) - V(old) 
V(old) 


< 0.03 


for all elements. (That is, all elements must satisfy this 
criterion to stop the iteration process.) A dummy subroutine 
(AITERY) is coded to allow for further modifications of to 

incorporate the utilization of perception thresholds in the 
display. 

The user must be careful to structure the matrix H , since 

c 

XCOV is the covariance matrix of Z and not Z. (If 

=1, the two matrices are equivalent,) 

(2) An arbitrary polynomial (up to order six) of the form 
N(s)/D(s) may be multiplied by any transfer function prior to 
obtaining a Bode plot. This provides additional plotting 
capability for the user. 

(3) The input switch FCNTRL=T (in the input, Group 1 

cards) causes the structure of to be slightly modified. 

Using this option, F^ is given by: 
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F = 
c 





u 

" pp - s 

0 


u 

1 

« 

1 

tc, 
tr 1 

Ph - KHh - Gh A 


u 

-Sa 


Pa 

11 


It can be shown that this structure effectively bypasses 
the Kalman filter and the operator now only contains the optimal 
controller equations. This is obvious by inspection by noting 


= [X^ x'^ X^], 


and above decouple$ the filter states X 


from the closed loop system. (The filter terms are still 
included, but are decoupled from the system.) This implementation 
was used to minimize changes to the structure of F^. 

2.3 INPUT SPECIFICATIONS 


The program receives its input through cards and disk files 
2.3.1 Card Input 


I ! 


There are two cards required for each run. Card 1 is a 
title card for output identification. Card 2 is the option 
specification card. All succeeding input is dependent upon the 
options selected on card 2. All matrices utilizing card input 
are read in by rows. The default format is 8E10.4. All matrices 
utilizing disk file input are read in by rows. The default 
format is 4E20.13. Table 2.2 describes both the required and 
optional card input. 

2.3.2 Disk Files 

Depending on which options the user elects, the program can 
use up to seven disk files (logical units 9-14, 21). 
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2. 3. 2.1 Reading F and G (Plant Dynamics) from Disk (Unit 9) 

The unit for reading F and G (plant dynamics) from disk is 
referred to as lUFG and has a default value of 9. This unit need 
only be assigned when DISK = T. 

lUFG is a formatted file. F and G must be written by rows. 
They are read by a variable format stored in array IFHT2(13). 

The default value is (4E20.13) and may be altered via NAMELIST/ 
OPTION/. • 

2.3. 2.2 Scratch Unit (Unit 10) 

Logical Unit 10 is a scratch file. It should always be 
assigned as a temporary file. 

2.3.2. 3 Output from Plant Dynamics (Unit 11) 

Output from plant dynamics calculations is stored on a file 
referred to as lUP and has a default value of 11. This file 
contains the characteristic equation and the numerator matrices. 
It should be assigned v/henever PLANT = T. It v/ill usually be a 
temporary file. If it is assigned as a permanent file, on a 
subsequent run, the user may set JENTP = 1 on Card 2, Group 2 to 
produce further transfer functions without redoing previous 
calculations. 

2. 3. 2. 4 Output from Human Operator (Unit 12) 

Output from the human operator is stored on a file referred 
to as lUH and has a default value of 12. It should be assigned 
whenever HUMAN = T. Its description and use is the same as for 
Unit 11, Section 2. 3. 2. 3. 
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(Note: If assigned as a permanent file, the user can use 

file lUH as if it contained output from plant dynamics, and use 
JENTP = 1 on succeeding runs to get additional transfer func- 
tions.) 

2. 3. 2. 5 Output from Human Operator + Plant Dynamics (Unit 13) 


Output from the human operator + plant dynamics is stored 
on a file referred to as IU3 and has a default value of 13. It 
should be assigned whenever HOPLT = T. Its description and use 
is the same as for Unit 11, Section 2. 3. 2. 3 (see note. Section 
2 . 3 . 2 . 4 ) . 

2. 3. 2. 6 Output from Closed Loop (Unit 14)‘ 

Output from the closed loop is stored on a file referred to 
as lUC and has a default value of 14. It should be assigned 
whenever CLOSED = T. Its description and use is the same as for 
Unit 11, Section 2. 3. 2. 3 (see note. Section 2. 3. 2. 4). 

2. 3. 2. 7 Writing PC, GC, HC, WBIG, VY, VA to Disk (Unit 21) 

The unit for writing FC, GC, HC, WBIG, VY and VA (if used) 
is referred to as lUSYS and has a default value of 21. It is a 
binary file and need only be assigned whenever SAVSYS = T. 



Table 2.2 
Card Input 


HU 

CARD 

NO. 

FORMAT 

' VARIABLE 
NAME 

TYPE 

DESCRIPTION 

1 

1 

13A4 

ITITLE 

A 

52 character title to be 
printed at the top of each out- 
put page 


2 

NAMELIST 

/OPTION/ 

PLANT 

L 

= T, required for plant dynamics 
calculations (Default = F) 




HUMAN 

L 

= T, required for human opera- 
tor calculations (Default = F) 




HOPLT 

L 

= T, required for human opera- 
tor and plant dynamics calcu- 
lations (Default = F) 




CLOSED 

L 

= T, required for closed loop 
calculations (Default = F) 




COVAR 

L 

= T, required to compute 
covariance (Default = F) 




MINOUT 

L 

= T, to obtain "minimal" out- 
put (Default = F) 




BACKWD 

L 

® T, to use backward Leverrier 
algorithm to compute adjoint 
matrix for no. of states < 5. 
(Automatically uses backward 
Leverrier algorithm if no. of 
states ^5). (Default = F) 



• 

DISK 

L 

= T, if plant dynamics F and G 
matrices are to be read from 
Unit 9. 

• 



SAVSYS 

L 

= T, to write FC, GC, HC, W, 
Vy, VA to Unit 21 in binary 
format (Default = F) 

1 

1 

1 ; " ' 

i ■ ! 

; : i 

i ’ ' 

i ' I 

■ j . i 


;JUNIT 

I 

Unit no, all matrices (except 
F and G for plant dynamics 
if DISK = T) and matrix 
elements will be read from 
(Default = 5) 




IFMTC13) 

A 

Array containing format to be 
used when reading matrices 
from JUNIT (Default = 

(8E10.4) 


Note: All Group and card numbers marked with * are conditional input. 
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Table 2.2 (Continued) 


GROUP 

NO. 

CARP 

NO. 

FORMAT 

VARIABLE 

NAME 

. 

TYPE 

1 

DESCRIPTION 

1 



IFMT2(13) 

A 

Array containing format to be 
used when reading F and G 
from Unit 9 


1 : 1 

!. 




(Default = (4E20.13)) 

i : 




ITERG 

L 

= T, for iterating on Vy 
and/or Va (Default = F) 

1 ' ■' 



ITERVY 

I 

No. of iterations to do on 
Vy (Default = 0) 




ITERVA 

I 

No. of iterations to do on 
Va (Default = 0) 




RH04 : 

D.P, 

Constant for Vy iteration 
(Default = .0314159 




RHOY 

D.P. 

Constant for Va iteration 
(Default = .092477) 

1 

1 


' 

FCNTRL 

L 

= T, to switch FC matrix to 
control only. (Default = F) 

■ 

2* 

■ 

1 

13A4 

ICMMNT 

A 

♦Include Group 2 only if 
PLANT = T 

Comment card for user identi- 
fication of start of plant 
dynamics input. 


2 

815 

NSP 

I 

No. of states (£30) 





I 

No. of controls (<4) 





I 

Np. of measurements (<_10) 




NOP 

I 

No. of noise (_<4) 




NTFP 

I 

No. of transfer functions 
desired (may be 0) 




JEIGFP 

I 

Eigenvalue flag for F 
=;-l, compute and skip all 
other plant dynamics calcu- 
lations 

= 0. do not comoute 

= 1, compute and continue 

calculations 




JFSNGP 

I 



Singularity flag for F 

j 


Note: All Group and card numbers marked with * are conditional input. 


24 











Table 2.2 (Continued) 


B 

CARD 

NO. 

— 

FORMAT 

VARIABLE 

NAME 

TYPE 

DESCRIPTIO^ 






= 0, unknown 
= 1 , F is singular 




JENTP 

I 

Entry flag 

= 0, regular run 

= 1, if adjoint matrices 

were previously computed and 

saved, and only specific 

transfer functions are 

desired 


3* 

615 



*Include NTFP sets of cards 
3-3B 




IRWCLP(l) 

I 

Row no. of transfer function 
des i red 




(2) 

I 

Column no. of transfer func- 
tion desired. If IRWCLP (2) 
<0, the transfer function 
produced will be (IRWCLP(l), 
(-IRWCLP(2)+1)*S. 




(3) 

I 

= 0, zeros only 
= 1, zeros, residues 
= 2, bode plots only 
= 3. all 




(4) 

I 

Frequency (in radians/second) 
of lower bound of bode plot 
(log) 




(5) 

I 

Upper bound of bode plot 
(log) 




(6) 

I 

= 0, normal size bode plot 
= 1, double size bode plot 
(magnitude scale only) 


3A* 

(I5,7E10.5) 

NCOEFN 

I 

No. of coefficients in the 
numerator of the polynomial 
(N/D) which multiplies the 
transfer function before 
plotting. (_<6) 




COEFNP 

D.P. 

Coefficients in descending 
powers of S 

♦Include only if IRWCLP(3)>2 


Note: All Group and card numbers marked with * are conditional input. 
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Table 2.2 (Continued) 


GROUP 

NO. 

CARD 

NO. 

FORMAT 

VARIABLE 

NAME 

TYPE 

DESCRIPTION 


3B* 

(I5,7E10.5) 

NCOEFD 

I 

No. of coefficients in the 
denominator of the polynomial 
(N/D) which multiplies the 
transfer function before 
plotting. (£6) 




COEFDP 

D.P. 

Coefficients in descending 
powers of S 

♦Include only if IRWCLP(3)^2 


4-4n* 

IFMT 

FP 

D.P. 

F matrix (NSP x NSP) 
♦Include only if JENTP = 0. 


5-5n* 

IFMT 

GP 

D.P. 

G matrix (NSP x NCP) 
♦Include only if JENTP = 0, 
JEIGFP f -1 


6-6n* 

IFMT 

HP 

D.P. 

H matrix (NMP x NSP) 
♦Include only if JENTP = 0. 
JEIGFP / -1 


7-7n* 

IFMT 

DP 

D.P. 

D matrix (NMP x MCP) 
♦Include only if JENTP = 0. 
JEIGFP f -1 


8-8n* 

IFMT 

GAMMP 

D.P. 

r matrix (NSP x N0P) 
♦Include only if JENTP = 0, 
JEIGFP f -1 

3* 





♦Include only if HUMAN = T 

i ' 

i i ! 

1 

13A4 

ICMMNT 

A 

Comment card for user identi- 
fication of human operator 
input 

1 ■ 

2 

715 

NSH 

1 

I 

No. of states (<15) 




kw 

I 

No. of controls (<_4) 




^MH 

I 

No. of measurements (<^10) 




NOH 

I 

No. of noise (£4) 




NTFH 

I 

No. of transfer functions 
desired (may be 0) 




JEIGFH 

I 

Eigenvalue flag (see Group 2 
card 2) 




JFSNGH 

I 

Singularity flag (see Group 2, 
card 2) 


Note: All Group and card numbers marked with * are conditional input. 
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Table 2.2 (Continued) 


GROUP 

NO. 

CARD 

NO. 

FORMAT 

VARIABLE 

NAME 

TYPE 

DESCRIPTION 

!\ ■ 

3* 

615 

IRWCLH(l-6) 

I 

See Group 2 , card 3 
*Include NTFH sets of cards 
3-3B 


3A* 

(I5,7E10.5) 

NCOEFN 

I 

See Group 2 , card 3A 



■ 

COEFNH 

D.P. 

See Group 2, card 3A 


3B* 

(I5,7E10.5) 

1 NCOEFD 

I 

See Group 2, card 3B 




1 COEFDH 

D.P. 

See Group 2, card 3B 


4-4n* 

IFMT 

'fh 

D.P. 

F matrix (MSH x NSH) 


5-5n* 

IFMT 

GH 

D.P. 

G matrix (NSH x NCH) 


6-6n* 

IFMT 

HH 

D.P. 

H matrix (NMH x NSA) 


7-7n* 

IFMT 

DH 

D.P. 

D matrix (NMH x NCH) 


8-8n 

IFMT 

GAMMH 

D.P. 

r matrix (NSH x N0H) 


9-9n 

IFMT 

QH 

D.P. 

Q matrix (NSH x NSH) 

4* 




i : 

♦Include only if HUMAN = T 


1 

13A4 

ICMMNT 

A' 

! 

Comment card for user iden- 
tification of augmented 
states input 


2 

15 

NSA 

I 

No. of augmented states 


3 

IFMT 

FA 

D.P. 

FA matrix (NSA x NSA) 


4 

IFMT 

GA 

D.P. 

GA matrix (NSA x NCH) 


5 

IFMT 

HA 

D.P. 

HA matrix (NCA x NSA) 


6 

IFMT 

DA 

D.P. 

DA matrix (NCH x NCH) 

'■ ■ 

f 

7* 

' 

IFMT 

VA 

* 

D.P. 

VA matrix (NCH x NCH) 

♦Include only if COVAR = T 

1 

I 5* 





♦Include only if HOPLT = T. 


1 

i 13A4 

ICMMNT 

A 

Comment card for user iden- 
tification of human operator 
+ plant dynamics input • 
section 


2 

15 

NTF3 




No. of transfer functions 
desired (may be 0) 


Note: All Group and card numbers marked with * are conditional input. 
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Table 2.2 (Continued) 


GROUP 

NO. 

CARO 

NO. 

FORMAT 

VARIABLE 

NAME 

TYPE 

DESCRIPTION 


3* 

615 

TRWCL3(l-6) 

j 

1 

I 

. 

See Group 2, card 3 
*Include NTFS sets of 
cards 3-3B 


3A* 

(I5,7E10.5) 

NCOEFN 

I 

See Group 2, card 3A 



' 

C0EFN3 

D.P. 

See Group 2, card 3A 


3B* 

(I5.7E10.5) 

NCOEFD 

I 

See Group 2, card 3B 




C0EFD3 

D.P. 

See Group 2, card 3B 

6 





*Include only if CLOSED = T 


1 

13A4 

i 

ICMMNT 

A 

Comment card for user iden- 
tification of closed loop 
input section. 


2 

hs 

NMC 

I 

No. of measurements (_<35) 




NTf C ! 

I 

No. of transfer functions 
desired (may be 0) 




JEIGFC 

I 

Eigenvalue flag 
See Group 2, card 2 




IDNTHC 

I 

= 0, HC matrix will be 
read in 

= 1 , HC matrix set to 
identify 


3* 

615 

IRWCLC(1~6) 

I 

See Group 2, card 3 
♦Include NTFC sets of card 
3-3B 


iA* 

(I5,7E10.5) 

NCOEFN 

5 I 

See Group 2, card 3A 




COEFNC * 

D.P. 

See Group 2, card 3A 


3B* 

(I5.7E10.5) 

NCOEFD 

I 

See Group 2, card 3B 




COEFDC 

D.P. 

See Group 2, card 3B 

i 

4 

IFMT 

CC 

D.P. 

Matrix for augmenting GC 
(NSA X NCH) 

! 

5 

IFMT 

BC 

D.P. 

Matrix for augmenting GC 
(NCH X NCH) 


Note: All Group and card numbers marked with * are conditional input. 



















Table 2.2 (Concluded) 


DESCRIPTION 


H matrix for closed loop 

(NMC X NSC, NSC = NSP + 

NSH + NSA). *Include only 
if IDNTHC = 0 

*Include only if HUMAN = T 
g matrix (NCH x NCH) 

W matrix (NOH x NOH) 

Vy matrix (NMH x NMH) 


: All Group and card numbers marked with * are conditional input 

= Alphanumeric 
= Integer 
= Logical 

= Double Precision 







2 . 4 PROGRAM FLOWCHART 

The following pages contain the flowchart for the transfer 
function progiam. 


i j 


ij 

y 


Li 

t 


? i 


I i 


i. i 
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O'-IGINAL PAGE IS 
•1? POOS QUALITY 
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COMPUTE CONTROL, 
FILTER GAINS. 


JEIGFH=0? 


COMPUTE 

EIGENVALUES 


COMPUTE 
TRANSFER FCN, 


NO YES 

— ^JEIGFH<0? >— 


WRITE ALL MATRICES, 
CHAR. POLY. TO FILE lUH 


COMPUTE ZEROES, 
RESIDUES, IF FLAGGED 


PRODUCE BODE PLOT, 
IF FLAGGED 


END B 
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JEIGFC=0? 



COMPUTE 

TRANSFER FUNCTION 


WRITE ALL 

MATRICES, CHAR. POLY. 
IP-FILEIUC 


COMPUTE ZEROES, 
RESIDUES IF 
FLAGGED 


COMPUTE BODE PLOTS 
IF FLAGGED 











2.5 TSS/360 OPERATION 


The following control cards are required to execute the TF 
program on the NASA ARC TSS/360 system. 

(1) L0G0N... 


(2) 

DDEF 

FTlOFOOl, 

VS, 

SCRATCH, DISP = NEW, RET 

= T 

(3) 

DDEF 

FTllFOOl j 

VS, 

PLANT, DISP = NEW, RET = 

T 

(4) 

DDEF 

FT12F001, 

VS, 

HUMAN, DISP = NEW, RET = 

T 

(5) 

DDEF 

FT13F001, 

VS, 

H0PLT, DISP = NEW, RET = 

T 

(6) 

DDEF 

FT14F001, 

VS, 

CL0SED, DISP = NEW, RET = 

^ T 

(7) 

DDEF 

FT21F001, 

VS, 

TIME. FI, DISP = NEW 


(8) 

JBLB 

NEWTFLIB 




(9) 

L0AD BL0CK$$ 





(10) CALL MJUN$$ 

(11) {data 

(12) L0G0FF 

The previous transfer function described in the HOCOM user’s 
manual (see Section I) is executed identically by the above 
sequence with the following exceptions: 

(a) eliminate card #8 

(b) replace card #9 by LOAD BLKDTl 

(c) additional input data required in #11 

i 

The description of the above execution deck follows: 

(1) Usual L0G0N card. 

(2) Scratch file used by transfer function. It is always 
needed and should always be assigned temporary 

(RET = T parameter) . 
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(3) “(6) Files needed depending on options chosen. I.e., if 
closed loop dynamics is not desired, card (6) can be 
eliminated. The RET = T parameter makes these files 
temporary. 

(7) If the closed loop system is to be saved for use by 
the time varying program, this file must be assigned. 

(8) Library of routines. 

(9) Load in block data. 

(10) Execute program. 

(11) Input cards. 

(12) Usual L0G0FF. 
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2.6 SAMPLE EXECUTION 



The following represents a sample execution (with descrip- 
tion and output listing). 



i t 

i I 
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1 

1 : r' . r I ' 

2 

' ' r 

3 

i • r-i' 

4 


5 

OL.'Sl'^ = 

6 

^ I d ~ 

7 

T -V v. -p 

8 

uE.jO 

9 

I . POP ?i''. 

10 

2 

11 

-2. 

12 

J . 

13 

1 

i • w 

14 

J . 

15 


16 


17 

J . 

18 

, p 

19 

i.'jp'jr F‘ji\ 

20 

4m 

21 

1 

22 

J 1 .i 

23 

! 1 1 .( 

24 


25 

1 . p 

26 

•; . 0 

27 

1 . p 

28 

4 . 0 

29 

0.0 

30 

1 . l; 

31 

1 

m V 

32 

o.c 

33 

0. a 

34 

AUG.'IENTEO 

35 

2 

36 

0 . i' 

37 

“266 . o~ 

38 

J . i> 

39 

44 . 444 

40 

la. 


- ■! r-LA. 1’ uY-jA IJCf 


J 


1 . 


1 


41 

42 

43 

44 

45 

46 

47 

48 

49 
50j 

51 

52 

53 

54 


J 


A’ 

J 

>1 


ICS 


a 


1 . 0E-:b 


1 . c 


J .r 


i . t’ 

■^ 26 .65 


J . d 

*•2. J 
3. 1 04d 

LOS ED LOOP IJPUr 


0 

1 
J 
J 


. J 
1 


1 

1 

{] 

0 


0 

3 


y. « 

J • (ij 

. 

AUOITTOiv'AL fiU-IAN INPUT 
C5 


a a 


123 


L PAGE IS 
UALITY 
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••• INPUT ••• 

ITITLE = SAMPLE hUK — TBAMSfEH EUNCTION 
THE OPTIONS HAVE BEEN SET AS FOLLOWS 
PLANT DXUAMICS :T 
HUHAN LKMAMICS = T 
CLOSEO LOOP KESPONSE = T 
UUi.AN ♦ PLANT UICNAMICS s F 
MIl.IliAL OUTPUT FLAG = F 

PACK'nARD LEVeRRIeB FOR NO. STATES LESS THAN 5 = 
ALL MATRICES WILL BE READ WITH FORMAT (8(E10.l») ) 
FRoli UNIT S 
COVARIANCE FLAG = T 
DISK HEAD FOR F AND C (PLANT) s F 
WRITE rC,GC,HC,WbIG,VI,VA TO DISK FLAG = T 
f AND G (PLANT) ON DISK USE FORMAT (J1(E20.13) ) 
UNIT FOR FC.GC,flC,NBIG,VI,VA = 21 
SWITCH To SET FC MATRIX TO CONTROL ONLY = F 
FLAG TO ITERATE OH VY AND/OR VA = F 
•••PLANT DYNAMICS INPUT 
NO. OF STATES =2 
MO. OF CONTROLS = 1 

NO. OF MEASUREMENTS = 1 

NO. OF NOISE SOURCES = 1 

EIGENVALUE FLAG FOR F = 0 

SINGULARITY FLAG FOR F = 0 

ENTRY FLAG =0 

0 TRANSFER FUNCTIONS WILL BE DONE. 

SAMPLE RUN TRANSFER FUNCTION 


F 


PLANT DYNAMICS ” 

1 2 

1 -2.0000D 00 0.0000 

2 1.0000D 00 l.OOOOD-05 


H 


PLANT DYNAMICS 
1 

1 0.0000 


2 

1.00000 00 


2 

O.OOCO 


GAMMA, PLT. DYN. INPUT 

1 

1 l.OOOOD 00 

2 0.0000 


D, PLANT DYNAMICS 
1 

1 0.0000 


G, AS INPUT 

1 

1 0.0000 

2 l.OOOOD 00 

SAMPLE RUN TRANSFER FUNCTION 


AUGMENTED G AND GAMMA MATRIX 



0.0000 
L..oooorj „.Cif^ 


1.0000D 00 




O 

o 


II. MAThlX CF TRAKSFch rU.MCTIO.NS 
•»»riUMAH DYNAiiXCS IHPUT 
I.O. Or bTATtS s ?. 

Or COUTliOLS = 1 

Or MEASU htiiilLhTS = 

Or (lOISt SObhCbS = 
EIuci.VALUfc FLAG FOh f : 
SliiOULAKITY FLAG FOh F 


THE FIRST 


1 COLUMNS = If OVER 0 


The second 


1 COLUMNS = 


i— 

Y OVER W 


nr:. 


I.U 

1.0 

C.U 


1 


1 ThAMSFEH FONCTIOMS MILL BE DOIfE. 
hO»* COL. OPTION LOG HIH LOG MAX PLOT FLAG 
t 1 3 - 1 1 0 


SAMPLE HUN TRANSFER FUNCTIO.N 



e, 

HUM 

DXH. INPUT 


0 



1 

2 



1 

-2.0000D 00 

0.0000 



2 

1.0000D 00 

1 .0000D-G5 

Q 

■ 

G. 

HUM 

DIN. INPUT 


0 


1 

1 

0.0000 




2 

1.OG0OD GO 







*-* 

H, 

HUM 

DYK. INPUT 


Q 



1 

2 


1 

0.00 3 

1.0000D 00 

d 





0. 

HUM 

. DYN. INPUT 
1 

0.0000 


© 


1 


o 

GAMMA , 

HUH. DYN. INPUT 
1 




1 

1.0000D 00 


Q 


2 

0.0000 


1 

Q 

Q, 

HUM 

. DYN, INPUT 




1 

2 


i 

1 

0.0000 

0.0000 

Q 


2 

0.0000 

1.0000D 00 

SAMPLE 

bun — TRANSFER 

FUNCTION 

O 

•••AUG 

■Ic-NTEO STATES INPUT. 

NO 

. OF 

STATES = 2 



NO 

. Or* 

CONTROLS = 1 


O 

NO 

• Or' 

MEASUREHENTS = 

1 


FA 


2 


0.0000 1.0000D 00 


MUM. COEF.t 10000E 01 
DEN. COEF.. 10000E 01 




o o 


tri Q 

o s 
o ^ 




1 




Ss 3 


i~ 


«• .. 


GA 


O.COOG 

CV 


1 2 
1.6000D 01 1.UOOOD CO 


- 2.00000 00 


1 1. 04800-01 

SAHPLt RUN ---_Tm«iSFER fUNCTION 

CLOSED LOOP RESPONSE INPUT 
NO. OF STATES = 6 

NO. OF CONTROLS =3 
NO. OF ilEASUREMEKTS = 7 

SLT HC - IDENTITI MATRIX FLAG = 1 

EIGENVALUE FLAG =0 

1 TRANSFER FUNCTIONS WILL BE DONE. 

ROW COL. OPTION LOG HIN LOG MAX PLOT FLAG 
113 -1 1 0 


NUM. COEF..10000E 01 
DEN. CbEF..10000E 01 


CC MATRIX 


0.0000 
0. 0000 


1.00000 00 


UNAUGMEHTED CLOSED LOOP HC 

1 2 3 ' 

1 1.0000D 00 0.0000 0.0000 ( 

2 0.0000 l.OOOOD 00 0.0000 ( 

3 0.0000 0.0000 l.OOOOD 00 ! 

4 O.OOQO 0.0000 0.0000 1.00( 

•j 0.0000 0.0000 0.0000 ! 

6 0.0000 0.0000 0.0000 ( 

••• end of INPUT. 

“i* PLANT DYMALICS OU'^PU;'' 

OA TlST Or The. iiU;-.EhIC ACCUhACV Ur Tnc. ..r. . i-.n !• x i:.i. 

SHOULD BE IDLNTICALLX ZERO ACCORDING TO Tclt ’^HiiOhl. 

XT WAS CALCULATED FDR THIS PhCjBLLd AS 


4 

0. OGGO 
0.0000 
0.0000 
,00000 00 
0.0000 
0.0000 


N-Tti B (iAThI 

1 ? 

•TED’'' ■ f ":"". r ■ 

2 -c.^‘il2DilT"’-6.^*JT3C-ir~' 

LE^'ERHIER algorithm was computed I.N’ THE FORWARD DIRECTION ONLJ. 


c o c> 




0.0000 
0.0000 
0 . 0000 
. 0000 
.0000 
OD 00 


I 

i 

I 





•ill* ■■ I ml «-« <n i« mi ,j*i m i 





tiu hfc« atO'k- (■ mT fi Isrr 


0 



— — — — + 

KUHEHATOR MATRIX COErFieiENT FOR S»» ?. 




c 


1 ?. 

\ 0.0000 o.oooc 

MUMERATOfc MATRIX COEt'FICIERT FOR S»« 1 


— --1 : ■ - • 2 . - 

1 1 .OQGOD 00 O.COOO 

RUHERATOR MATRIX COEFFICIENT FOR S»* 0 


1 2 

1 2.0000D 00 1.00000 00 

MATRIX OF DC GAINS 


1 2 

1 -1.00000 05 -5.000CD 04 

SAMPLE RON — TRANSFER FUNCTION 


' POLES OF THE TRANSFER FUNCTION 
+ 

COEFFICIENTS Or THE CHARACTERISTIC EQUATION IN DESCENDING POWERS OF S = 

1.00000D 00 1.99999D 00 -2.00000D-05 

•POLES OF THE TRANSFER FUNCTION AS COMPUTED FROM THE ABOVE POLXNOMIAL = 

0 (-2.00000D 00) + J( 0.00000) 

HAG, ANGLE, ZETA = 0.200GD 01 0.1600D 03 0,|000D 01 

1 ( 1.00000D-05) + J( O.OCOOO) 

KAG, ANGLE, ZETA = 0.1000D-04 0.0000 -0.1000D 01 • 

SAMPLE RUN TRANSFER FUNCTION 

•••» CASE NO. 1 


GSML 

1 

1 

5.0000D-02 

WBIG 

1 

1 

1 . 7000D 01 

Vt 

1 

1 

1 .2300D-03 


SAMPLE RUN - — TRANSFER FUNCTION 
CONTROL AND FILTER GAINS. 



P 


§1 


1 ■ ? 

1 6.90960-01 1). 47210 00 


FILTER GAIN 

1 

1 9.06^60 01 

2 1.34640 01 

SAKPLt fiUH — TfiAnSFER FUNCTION 

••• HUMAN OPERATOR 07NAHICS OUTPU1 
MATRICES FOR HUMAN DYNAMICS. 


FHH 


•2.0000D 00 -9.06360 01 
2.36200 00 -4.61940 00 


0.0000 


0.0000 


3 

0.0000 
1.SOOOO 01 
0.0000 


4 

0.0000 
1.60000 00 
1 .OOOOD 00 


•3.0710D 01 -1.9U76D 02 -2.66670 02 -2.66670 01 


GHH 


9.06360 01 
1.34640 01 
0.0000 
, 0.0000 




HHH 

12 3 4 

1 1.36200 00 6.94430 00 1.60000 01 1.60900 00 

A TEST OF THE NUMERIC ACGURACl OF THE ONE-WAY LEVERRXER METHOD IS THAT THE N-TH B MATRIX 
SHOULC BE IDENTICALLY ZERO ACCORDING TO THE THEORY, 

IT WAS CALCULATED FOR THIS PROBLEM AS — - 


N-TH B MATH I 


1 


1 .1.36420-11 

2 -6.4570D-13 

3 2.7265D-12 

4 -9.0949D-16 

LEVERRIEft ALGORITHM WAS COMPUTED 
SAMPLE RUN TRANSFER FUNCTION 


2 

-2.9104D-11 
-3. 6360D-11 
- 3 . 36610-11 
6.46130-10 


1 .7621D-10 
3.4661D-1 1 
3.6360D-1 1 
1.4S62D-1 1 


-9.6634D-13 
2. 0916D-12 
9.0949D-13 
-2. 4S56D-1 1 


IN "^KE FORWARD DIRECTION ONLY. 


NUMERATOR MATRIX POLYNOMIAL OF THE TRANSFER FUNCTION 


NUMERATOR MATRIX COEFFICIENT FOR S*» 4 


1 0.0000 

iU.-i£;hATOfi MATRIX COEFFICIEI''' FOR E*‘ 3 


1 2.4566D 02 

her M A COF'’ =’’--CIEI"^ -R-Ofl -3 




iTH” 




i.UHEjiA'tUri MAfhiX’cOifrri^IEijV tiih S’-» i’ 






. 1 -'4.677SU 04 

;*Ui‘iti.AT0R MATRIX COErFICIENT FOR S*» 0 


1 -1.40190 05 

MATRIX OF DC GAINS 


1 -B.2^95D 00 

SAMPLE RUN THAMSFER FUNCTION 


■ POLES OF THE TRANSFER FUNCTION 

COEFFICIENTS OF THE CHARACTERISTIC EQUATION JN DESCENDING POWERS OF S = 


l.OOOOOD 00 3.31B6RD 01 1.02S22D 03 7.01979D 03 1.70357D 04 

POLES OF THE TRANSFER FUNCTION AS COMPUTED FROM THE ABOVE POLINOMIAL = 


cn 


0 (-1 

.25319D 

01 ) 

+ J( 2.53228D 

01) 




MAG. ANGLE, 

ZETA = 


0.2825D 02 

0. 1163D 

03 

0.9933D 

00 

1 (-1 

.25319D 

01) 

+ J(-2.53226D 

01 ) 




HAG. ANGLE, 

ZETA = 


0.2825D 02 

0 . 11630 

03 

0. 94350 

00 

2 (-9 

.06lblD 

00 ) 

+ J( 2.20050D 

00) 




HAG. ANGLE, 

ZETA = 


0.9620D 01 

0. 1516D 

03 

0.8793D 

00 

3 (-9 

.06161D 

00) 

+ J(-2.2005CD 

00) 




MAG, ANGLE, 

ZETA = 


0.9620D 01 

0» 1516D 

03 

0.6793D 

00 


SAMPLE RUN TRANSFER FUNCTION 

♦ + 

' ZEROES OF THE TRANSFER FUNCTION 

ZEROES OF THE TRANSFER FUNCTION BETWEEN Y( 1 ) AND 0( 1) = 

1 ( 1.99992D 01) + J( 0.00000) 

HAG. ANGLE, ZETA = 0.2000D 02 0.0000 -0.1000D 01 

2 (-6.66657D 00) + J( 0.00000) 

MAG, ANGLE, ZETA = 0.6667D 01 0.1BOOD 03 0.1000D 01 

3 (-'(.2O005D 00) ♦ J( 0.00000) 

MAG, ANGLE, ZETA = 0.il280D 01 0.1E00D 03 0.1000D 01 

■COEFFICIENTS 0? THE NUMERATOR POLXNOHIAL IN ASCENDING POWERS OF S = 


o o 

Eg pj 
Pi 

i-d p 

§s 

§ ^ 


p 


K c3 


I 


-1.^01950 05 -il.RTTSOD 04 -2.224030 03 

RfcSlDUES Or THE TRANSFER rU.NCTI0.'i = 

POLE LOCATION 


2.«56aOD 02 


RESIDUb AT •’■HE POLE 


COMBINED fi£SIDUES=llAG»CGStwD»T + PHI) 
MAGNITUDE Pril(OEG) 


(-1.25314D 01) ♦ J( 2.53228D 01) ( 1.35698D 02) + J( 1.«036lD 02) 

(-4.06181D 00) J( 2.20050D 00) (- 1 . 285b« D 0 1 ) + J(-6.3‘>190D 00) 

FUNCTION TO BE PLOTTED IS I ( 1) TO U ( 1 ) 
rriARSFER function to be plotted 


3. 909630 02 
2.66635D 01 


-45.96773 

153.71120 


r n U ' p c 7 r T K w ^ 

f . 


r M M !•' 'ir nr j.J c, .fiicriER 


• V 


.,*T — y,/^/‘/'sur'C — 1. v>-*/ 

L£i'.'u;ilr. ATO It CUt Kr ICItl. T3 It. Li bSC :i(j 1 1 iiu ulii/Lh 
O.lOCGuOOt: 01 0.3^1o6AU OP O.IOPSPPPi'. O't G.7019V9P.-. G'l 0. 
LOGi. r.LOT vALUto - — 

Fi-ECu-./Ci .lACii.II'UDE PliASt; 

0.1OOQOE 00 0.1b309t. OP 0.1? ‘>‘>a O'il 
0.10<i71E 00 0.16i09t: 02 0.179*3^11 0^- 
0.10965E 00 0.18309E 02 0.179*3lE 03 
O.llUbPii 00 0.18309E 02 0.179't8E 03 
0.12G23E 00 0.18309E 02 0.179‘I6E 03 
0,. 125898 00 0.18309E 02 0.179'(3E 03 
0..13183E 00 0. 16310E 02 O.Wg'tlE 03 
0* 1300<1E 00 0.18310E 02 0.17938E 03 
0.1845‘IE 00 0.18310E 02 0.17935E 03 
0»15156E 00 0.18310E 02 0.17932E 03 
0,15o49b 00 0.1b311E 02 0.17929E 03 
0.1t596£ 00 0.18311E 02 0.17925E 03 
0.17378E 00 0.16311E OP 0.17922E 03 
0.1d197E 00 0.18312E 02 0.17918E 03 
0.19055E 00 0.1U312E 02 0.1791*18 03 
0.19953E 00 0.18313E 02 0.17910E 03 
0. 208938 00 0.18313E 02 0.17906E 03 
0.21878E 00 0. 1831*1E 02 0.17901E 03 
0.22509b 00 0.18314E 02 0.17897E 03 
0.239bfeb 00 0.18315E 02 0.17892E 03 
0.25119E 00 0.18316E 02 0.17887E 03 
0.26303E 00 0.18316S 02 0.17881E 03 
0.2754211 00 0.18317E 02 0.17876E 03 
0.2fc840E 00 0.18318E 02 0.17670E 03 
O.302GOE 00 0.18319b 02 0.17864b 03 
0.31623b 00 0.18320E 02 0.17857E 03 
0.33113E GO 0.18322b 02 0.17850E 03 
G. 34674E 00 0. 18323E 02 0.17b43E 03 
C*' 0.363G8E 00 0.1832*)£ 02 0.17636b 03 

0.3O019E 00 0.18326b' 02 0.17828E 03 
0.39b11E 00 0.16326b 02 0.17819b 03 
0.41687E 00 0.18330E 02 0.17811E 03 
0.43652E 00 0.18332E 02 0.17602E 03 
0.A57O9E GO 0.1833*18' 02 0.17792E 03 
0.47b63E 00 0.1b337E 02 0.17782b 03 
0.50119b 00 0.18339b 02 0.17771E 03 
0.524b1E 00 0.183*128 02 0.17760E 03 
0.5495*tE 00 0.183*168 02 0.177*198 03 
0.57544E 00 0. 183*19E 02 0.17736b 03 
0.60256b CO 0.18353E 02 0.17723b 03 
0.630968 00 0..18357fc 02 0.17710b 03 
0.6t)069b 00 0.18362b 02 0.17t95k. 0 3 
0.091838 00 0. 163678 02 0. 176808 03 
0.724448 00 0.16373b 02 0.17664E 03 
0.75858E 00 0.18378E 02 0.176*188 03 
0.79*1338 GO 0.183858 02 0.176308 03 
0.83176E GO 0.18352b 02 0.176118 03 
O.C70568 00 0.18400c 02 0.175918 03 
0.912018 00 O.IH'iObb 0? 0.17570b 03 
0.95459b 00 0.184178 02 0.175488 03 
0.100008 01 0.184268 02 0.175258 03 
0.10'<718 01 0.184378 02 0.175008 03 
0.109658 01 0.184488 02 0.17474b 03 
0.114828 01 0.184608 02 0.174468 03 
0.120238 01 0.184738 02 0.174168 03 
0.125858 01 0.184878 0? 0.173858 03 

5.1'j' , OV -'‘850: ‘ 2 0 • 5 1b • ...... , 

O.nTo'ib 0rT.l8517b u? 0.1/3158 03 . . : - 

0.14454c, 01 0. 185348 02 0. 17277b 03 

\ w- ' 


7G‘5";758 G5 



i:: 

ii' 


0 . 1*50^9^ 

^ ^ I 

0. Ic197t, 
0. 190951:. 
0 . 19 9*’ Tt 

0.20fc9^ii 
O.?l67ot 
0.22909c, 
0 . 2 >9 9^ c. 
0. 291 19c, 
0 . 2 6 1 0 H ii 
0 . 27992E 
0.2fio90t 
0.30199b' 
0.3U23E 
0.33113!:. 
0. 39674E 
0.3630fafc, 
0. 3b019t 
0. 39c1 ie 
0.416b7H 
0.43C92t. 
0.95709E 
0.47«(>3E 
0.901 19b 
0.929B1E 
0.999‘>4E 
0.97*i99B 
0.6029 6E 
0. 63096 b 
0 . o6069b 
0.69163c. 

0 . 72'49 36 
0.79b9bE 
0.79433c. 
0.03176E 
0.670966 
0.91201E 
0.99999b 

0 . loooob 

bOCc. PLOT . 


01 0.16970L 0? 


01 

01 

01 

01 

01 


01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
G 1 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
0l 
01 
01 


i 6 6094j»’"sit 


lot) 3oc. 02 
1B692b 0? 
Ib67*tb 02 
1b697E 02 
01 0.1b719b 02 

01 0.16742D 02 
0.1b76<lb 02 
0.1b769L‘ 02 
O.lboObb 02 
0.1bb29E 02 
b.lbdll2E 02 
0.l6fa97E 02 

1b869E 02 
1bb79E 02 
1b8faHE 02 
1bbb6E 02 
IbbbUE 0? 
l8b7bE 02 
l6b67E 02 
lbb92b 02 
lbb33b 02 
IbfalOE 02 
0.187.36E 02 
0.1b799E 02 
0.18724E 02 
lfa693E 02 
18661E 02 
18631b 02 
18603b 02 
0.1897EE 02 
0.16998E 02 
0.1 894 ME 02 
0.1b937E 02 
0.18939b 02 
0.18990E 02 
O.- 18 97 1b 02 
0.1860 ME 02 

02 0.18690E 02 
..MAGHITUOb OF 


0. 17194E 
0 . ^ b 

0 .^4'*f»K'99E 
0. 170M7E 
0. 16991b 
0. 16932E 
0.1 6869 E 
16801E 
16730b 
16693c. 
16973E 
16Mo7E 
16396b 
16300b 
16199E 
16093E 
0. 19982E 
0. 19869E 
0. 197Miib 
0 . 1 9 6 1 8b 
0. 19M88E 
0. 1939ME 
0. 19216E 
0. 1S079E 
0. 1M931E 
0. 11»78ME 
0. 1M639E 
0. 1U483E 
0 . 1 4 330E 
0,1817'IE 
0.1 MO 17 b 
0. 1389BE 
0. 13696E 
0, 13932b 
0. 13366E 
0. 13196E 
0. 13022E 
0. 128*tJlE 

if. 12660E 

0. 12471E 
0. 12279b 
KESP0.3SE 


03 


u 1 
03 
03 
03 
03 
03 
03 
0 3 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 
03 


t •- 
\ — 


- 1 - ■ 

u 


t:. 


\ 


20 . 




2. E~01 


1 .E-01 

bODc PLOT. . . PHASE Of TtlE BE3P0;'»S£. 


9. E-01 1 .E 00 

KKEui; EliCV , !. ( iiA O/SbC ) 


2. E 


00 


T.iE 00 


1 .E 01 


180, 


«i « < 


O'^IGINAL PAGE IS 
■ m QUALITY 




Ln;ii'tnhlb.l( ALCOhlTriK «AS COMHUTtiD IlJ '^Hb. !•’ 0 ii w A li D 

c 

- i * • V .■ t t- ■ I i,....: ; . 

+ — — 

' iNUhtfiATOH HAThxX fGLVNOHIAL 0£-rH-ii ’•ftANSFi-Ji. I 

+ — — — — 

N U M t. ti ATO H MATnIX COtKf'ICIcwT FuH S** 6 



1 

2 

3 

1 

0.0000 

0 . 0000 

0.0000 

2 

0.0000 

0.0000 

o.'oooo 

3 

0.0000 

0.0000 

O.OGOO 

6 

0.0000 

O.OQOO 

0. 0000 

5 

0.0000 

0 . 0000 

0.0000 

6 

0.0000 

0.0000 

0 . 0000 

7 

0.0000 

0.0000 

0.0000 


NUHgHATOfi MAThlX COEFFICIENT FOR S»» 5 


1 2 3 

1 0.0000 2. <15660 02 6.96630 00 

2 1.00C0D 00 0.0000 O.OGOO 

3 0.0000 O.OGOO 1.00000 00 

6 0.0000 9.0636D 01 0.0000 

5 0.0000 1,36660 01 1.00000 00 

6 0.0000 0.0000 0.0000 

7 0.0000 0.0000 0.0000 

NUMEhATOK MATRIX COEFFICIENT FOR S»» 6 


1 2 3 

1 0.0000 -1.7327U 03 3.6950D 01 

2 3.31660 01 0.0000 0.0000 

3 1.00000 00 2. *<5680 02 6.6131D 01 

6 0.0000 1.7676D 03 0.0000 

5 0.0000 6.28780 02 6.61310 01 

6 0.0000 0.0000 0.0000 

7 0.0000 -5.65950 03 -1.9876D 02 

NU.'icliATOh MATRIX CCEFFICIEHT FOR S»* 3 
SAMPLE RUN — - TRANSFER FUNCTION 


12 3 

1 2.65681) 02 -5.12230 06 -1.96690 03 

2 7.77*160 02 0.0000 0.0000 

3 3.31861) 01 -1.73270 03 8.8087D 02 

6 9.06360 01 3.61910 06 1.13690-13 

5 1.3666D 01 6.25920 03 8.8067D 02 

6 0.0000 -5.65950 03 -1.98760 02 

7 0.0000 -3.62860 06 -3.67110 03 
NUilEKATOR MATRIX COEFFICIENT FOR S»» 2 


1 2 3 

1 -2.22601) 03 -2.3376D 05 -3-6630D 06 

2 9.2638D 03 0.0000 0.0000 

3 1.0232U 03 -5.12230 06 8.8520D 03 

6 1.6063D 03 1.78060 05 -7.2760D-12 

■V 



6 0.0000 ~^.k?.obb 0« -1. '171 10 0 1 

7 -‘j. 459‘ib O'l O'! -2. 95 HU 0" 


NU.Sc.hATOfi I'lAThlX COc.KFICIE.4r FOfi S*» 0 


2 5 

2.b039D 00 -2.b039U 05 
0.0000 0.0000 
2.60390 05 6.27010 04 

2.16170 00 -7.27600-11 
2.93010-01 6.27010 04 

4. 67340-01 -4.6734D 04 
0.0000 0.0000 


1 2 3 

1 -5.00000-01 1.00000-05 -1.00000 00 

2 5.0000D-01 0.0000 O.OOOO 

3 6.07570-02 -1.00000 00 2.23620-01 

4 3.65460-01 -7. 70950-06 -2.59490-16' 

5 5.22510-02 -1. 04500-06 2.23620-01 

6 -6.33370-02 1.66670-06 -1.66670-01 

7 0.0000 0.0000 0.0000 

SAMPLE RUti TKAN.SFER FUNCTION 

NU.'-iEhATOR MATRIX COEFFICIENT FOR S«» 1 


1 2 3 

1 -4.67750 04 -2.80390 05 -2.05121) 05 

2 6.J18J10 04 0.0000 0.0000 

3 7.01980 03 -2.33740 05 4.58660 04 

4 3.49780 04 2.1617D 05 -4.36560-11 

5 5.05540 03 2.93010 04 4.58680 04 

6 -5.45950 03 -4.67330 04 -2.95140 04 

7 -2.33670 04 4.67340-01 -4.67340 04 

SAMPLE HUM . TRANSFER FUNCTION 

+----- — — + 


POLES OF THE TRANSFER FUNCTItuN 


COEr r rCIE.NTS OF THE CHARACTERISTIC EQUATIOIl IN DESCE.NDING POWERS OF S = 

1.000000 GO 3.518640 01 8.439150 02 1.079890 04 8.229830 04 2.678160 05 2.603690 05 

POLES OF THE TRANSFER FUNCTION AS COMPUTED FROM THE ABOVE POLYNOMIAL = 

0 {-7.0B451D 00) ♦ J( 1.689310 01) 

MAG, ANGLE, ZETA = 0.1832D 02 0.11280 03 

1 (-7.064510 00) + J(-1. 669310 01) 

MAG, AiJGLt, ZETA = 0.16320 02 -0.1 1280 03 

2 (-7.73166D 00) + J( 7.601430 GO) 

MAO, ANuLe, ZETA = 0.10640 02 0.13550 03 

3 (-7.731660 00) + J(-7.601'i3:J 00) 

HAG, ANGLE, ZE'^A = 0.1064U 02 -0.13550 03 

4 (-3.553690 00) ♦ J( 0.00000) 

MAG, ANGLE, ZETA = 0.35540 01 0.1800D 03 

5 (-2.00GOOD 00) ♦ J( 0.00000) 

MAO. ANCLE, ZETA - 0.20000 01 0.18000 03 

'•ItiPl}’’ FUC'-’-' ON , .... 


0.38670 00 
0. 36670 00 
0.71310 00 
0.71310 00 
0 . 10000 01 
0.10000 01 



1 

1 -1.40190 05 

2 1.4019D 05 

3 1.70360 04 - 

4 l.OfiOfaO 05 - 

5 1 . 465 ID 04 - 

6 -2.33670 04 

7 0.0000 
MATRIX Of DC GAINS 





'kotg 3 The . ^i.sr^ - 3'UkC'* ■ ‘ btf -li r ' - Aijr— .. j-, d ^ - . .._. 

1 ( 1.9909J!D 01) ♦ J( 0.00000) 

MAG, AtiGLt. ZtTA = O.POOOU OP 0.0000 -0.10000 01 

P (-6.66667D 00) + J( 0.00000) 

fiAG, AriuLt, Zf^A = P.6067& 01 O.TfcOOD OS 0.10001; 01 

S (-M.PoOGSt/ 0 0) * J( 0.00000) 

MAG, A«GLE, ZbTA = 0. iJPdOD 01 0.1800D OS O.IGOO!) 01 

COc-cflCItOTS OF THfc jlUl-iEKATOh PaLVMOMIAL II. ASCii.JUIlJG POWEHS OF S = 


-l.ltOlQ'iD GS -A.677‘iOG 04 -P.??40S0 OS 

htSlOUES OF Tilt TnAIJSFEfi rUHCTlOfU 

POLE LOCATION f 


2.4Se«0D 02 


HESIDUE AT THE POLE 


COMBINED RESIDUES=MAG*COS(kD»T+PHI) 
MAGNITUDE PHr(DEC) 


(-7.034‘>1D 00) ♦ J( 1.689S1D 01) 
{-7.731b6D 00) + J( 7.601430 00)' 

-3. *5^3690 00 

-2.00000D 00 

FUNCTION TO LE PLOTTED IS I( 1) TO U( 1) 
THAKSFEB FUNCTION TO BE PLOTTED 


(-6.734B6O-01 ) ♦ J(-6.933‘)30-01 ) 1.93321D 00 

( 1.14165D 00) ♦ J( 1.09210D 00) 3.1S978D 00 

3.757790-01 
-1.3121PD 00 


134.16723 

-43.72915 


NUMEBATOfl COEFFICIENTS IN DESCENDING ORDER 
G.P456795E 0 3-0 . 2224 027E 04-0 . 4677496E 05-0 . 1401 949E 06 
DEKOHIMATOH COEFFICIENTS IN DESCENDING ORDER 

0.1000000E 01 0.3516642E 02 0.8439154E 03 0. 1079889E 05 0.6229B26E 05 0.P67616OE 06 0.2B03895E 06 
BODE PLOT VALUES 

FhcOUENCi' MAGNITUDE PHASE - 


. 1 0 0 0 0 E 
. 10471E 
. 10965E 
. 1 1462E 
. 12G23E 
. 12509E 
. 1 3 1 6 3 E 
. 1.3B04E 
. 14454E 
. 15136E 
. 15649c 
. 16596E 
. 1737bE 
. 18197E 
. 1 905 5£ 
. 19953E 
.206936 
•21676E 
•22909E 
.23968E 
•25119E 
•26303E 
. 275UPE 
, 2 b o 4 0 c 
. 30200c 
, 3 162 3E 
.33113b 
. 3 4 6 7 4 E 
. 3 6 3 G 6 c 
, 36019E 
. 39t11L 
, 4 16&7E 
s b P 1. 


I.60312E 01 
I. 60323E 01 
I.60334E 01 
I.60346E 01 
I.60360E 01 
I. 60375E 01 
i. 6 0 39 IE 01 
I.60409E 01 
i. 60428 E 01 
I.60450E 01 
I.60473E 01 
.60499E 01 
.60527b 01 
.60558E 01 
.60591E 01 
.60626E 01 
.60669E 01 
.607136 01 
.60762b 01 
.60615E 01 
.60874E 01 
.60937E 01 
.610G7E 01 
.61064b 01 
:61166b 01 
.612596 01 
.613596 01 
.61469E 01 
.615696 01 
■61719E 01 
.61662c 01 

.62019b 01 
APlr.t, K (i1 


0.17644E 03 
0. 17627E 03 
0. 17610E 03 
0.17592E 03 
0.17572b 03 
0,17552c 03 
0.17531E 03 
0. 17509E 03 
0.17466E 03 
0.17462c 03 
G.17437E 03 
0. 174 ICE 03 
0.17362E 03 
0.17353c 03 
0.17323E 03 
0.17291E 03 
0. 17258c 03 
0.17224c 03 
0.17187E 03 
0.17149b 03 
0.17109E 03 
0.17 0666 03 
0.17 0 2 46 03 
0.16979b 03 
0 . 1 69 3 16. 0 3 
0. 16tr,2c 03 
0.16830E 03 
0.167756 03 
0.16719b 03 
0.16659E 03 
0. 165 9 6b 03 
0.16533E 03 
(5 1 0 3 


o 

O B 
O 

ciy 

C3 






Ol 

r-o 


0.7^^'i'it: 00 
0,7^Bf>0t; 00 
0.794^11; oc 
O.U3i?6k: 00 
0:&7096li 00 
O.9V?01b- oc 
0.9*>>i99b 00 
0.10000b 01 
0.10<17lb 01 
0. 10905b 01 
0.114U2L- 01 
0.1202^b 01 
0.12509b 01 

O.ni03b 01 

0. 13004b 01 
0.14454b 01 
0.15130b 01 
0. 15049b 01 
0.16590C 01 
0.17370b 01 
0.10197b 01 
0.19055b 01 
0.19953b 01 
0.20093b 01 
0.21070b 01 
0.22909b 01 
0,23900b 01 
0.25119b 01 
0.26303b 01 
0.27542b 01 
O.2004Db 01 
0.30199E 01 
0.31623b 01 
0.33113b 01 
0.34674b 01 
0.36300b 01 
0.30019b 01 
0.39011b 01 
0. 41667b 01 
0.43652b 01 
0.45709b 01 
0. 47063b 01 
0.50119b 01 
0.52401E 01 
0.54954b 01 
0.57544b 01 
0.60256b 01 
0.63096b 01 
0.66069b 01 
0. 69103b 01 
0.72443b 01 
0.75050b 01 
0.79433b 01 
0.03176b 01 
0.07096b Cl 
0.91201b 01 
0.95499a Ol 
0. 10000b 02 


0. 65447b 01 
0 . 6591 ^ 11 :: 01 
0.66420b 01 

0.66979b 01 
0.67575b 01 
O.602l0b 01 
0.60912b 01 
0.69650b 01 
0.70459b 01 
0.71320b 01 
0.72242b 01 
0.73227b 01 
0.74279b 01 
0.75400b 01 
0,76591b 01 
0-77055E 01 
0.79192b 01 
0.00604b 01 
0.02092b 01 
0.03655b 01 
0.052.93E 01 
0. 07004b 01 
0.00706b 01 
0.90641b 01 
0,92561E 01 
0.94544b 01 
0.96507b 01 
0.90603b 01 
0.10003b 02 
-0.10301b 02 
-0. 10524E 02 
-0.10749b 02 
-0. 10976b 02 
-0.11205E 02 
-0. 114 35b 02 
-0.11664b 02 
-0. 1 1892b 02 
-0.12119E 02 
-0. 12343b 02 
-0. 12565b 02 
-0.12783b 02 
-0. 12997E 02 
-0.13207b 02 
-0.13412b 02 
-0. 13613b 02 

-0. naiob 02 

-0.14003b 02 
-0. 14193b 02 
-0.14380b 02 
-0. 14566E 02 
-0. 14752b 02 
-0. 14939b 02 
-0. 15129b 02 
-0.15324b 02 
-0.15526E 02 
-0 . 15736b 02 
-0. 15956b 02 
-0, 16168b 02 


0. 15508b 03 
0.15399b 03 
0.15265b 03 
0. 15 168b 03 
0 . 15046b 03 
0.14920b 03 
0. 14790b 03 
0.14656b 03 
0. 14517b 03 
0.14374b 03 
0.14227b 03 
0.14075b 03 
0.13920b 03 
0.13761b 03 
0.135'97b 03 
0.13430b 03 
0.13259E 03 
0.13085b 03 
0.12907b 03 
0.12726b 03 
0.12542E 03 
0.12354b 03 
0.12164b 03 
0.11972b 03 
0.117/6E 03 
0.11578b 03 
O.inTOE 03 
0.11175b 03 
0.10970E 03 
0.10762b 03 
0. 10552b 03 
0. 10338E 03 
0.10121b 03 
0.99012b 02 
0.96770E 02 
0.94484b 02 
0.92146b 02 
0.89751E 02 
0.87291E 02 
0.64756E 02 
0.82137E 02 
0, 79423b 02 
0.76604b 02 
0.73666E 02 
0.70598b 02 
0.673b5E 02 
0.64015b 02 
0.60472b 02 
0.56745b 02 
0. 52818b 02 
0.4866ab 02 
0.44318E 02 
0.39721b 02 
0.34679b 02 
0.29784b 02 
0.24428b 02 
0. 18802b 02 
0.12900b 02 


• CODE PLOT. . .MAONITUOb Ob libSfOHbb 






f »•< 



, ‘It . I , I -.. . I i'l nn', ui ■ . I o(ii<ru 'j/ 

O.lijOOOt Ot>-0. Ulbaii 0? 0.1?900E 02 
liOOli PLOT . . .ilAUNITUUli OK HKSI'ONSK 



dz 


c 


u . + 


_20. ♦ 

l.K-01 

2.E-01 

*».E-01 

1.E 00 

2.E 00 

S.t 00 

1. 

boot HLOT . . . I'HASE 

ItO. •»»•«» 

OK THE fibSPOliSE 

• «_ ♦ 

KhEOUEIICY , 

, W (RAO/SEC) 





Cn 



0. 



•i.K-OI 1.K 00 ?.t 00 

PliEOUKIICY, h (HAD/SEC) 


S.E 00 


Tilh O’^I-.AOY COVnilIAIiCE HATiaX IS 


C«! ♦ 


1 .t-01 


?, li-01 


?.E 00 


00 


l.E 01 




s.e-01 1.E 00 

FIIEOUEHCY, b (HAU/SEC) 


THE :rrhAUt COVAKTANCE MAThlX IS 


E'-' 


1 P 3 H S 6 7 

i o.067bU 00 OO 00 -4.83950-01 8.18860-01 ?.?1670 00 

:> -3.?t»;^9U 00 4.;;5COO 00 4.83550-01 ?,5?61D 00 3.7?070-01 -4.67670-01 -9.35350-01 

3 -4.1)3550-01 4.83550-01 8.65470-02 3.72070-01 6.99860-02 -5.50870-02 -3.51190-01 

4 -3.28290 00 2..526I0 00 3.72070-01 2.52610 00 3.72070-01 -4.67670-01 -9.35350-01 

5 -4.03550-01 3.72070-01 6.99860-02 3.72070-01 6.99860-02 -5. 50970-0? -3.51190-01 

0 8. 18860-01 -4.67670-01 -5.50870-02 -4.67670-01 -5.50870-02 1.08770-01 -1.63760-15 

7 2.21670 00 -9.35350-01 -3.51190-01 -9.35350-01 -3.51190-01 -1.55430-15 3.69470 00 


COST IHOEX 

ItJX - -0.328290 01 

itJll = -0.327400 01 

liJ : -0.855691) 01 

•»» KC, GC, ilC, IU*IG, vy, ANI) VA KlilT-^EM TO DISK. 







Table 2.3 
Description of Sample Execution 


OBIGINAL PAGE IS 
OF POOR QUALITY- 


LINE 

NUMBER(S) 

GROUP 

TYPE 

CARD 

TYPE 

DESCRIPTION 

1 

1 

■■ 

Title for run 

2-8 



Option card. Plant, human and closed loop 
dynamics will be done. Covariance will be 
calculated and closed loop system will be 
$aved on file for use by time varying program. 

9 

2 

1 

Comment card indicating start of plant 1 

dynamics input 

10 


2 

Indicates 2 states, 1 control, 1 measurement, 

1 1 noise, no transfer functions, no eigenvalues, 
singularity of FP unknown, regular run 

11-12 


4 

FP matrix 

13-14 


i 5 

GP matrix 

15 


,6 

HP matrix 

16 


7 

DP matrix 

17-18 


8 

TP matrix 

19 

! 

3 ‘ 

1 

Comment card indicating start of human dynamics 
input 

20 

i 

! 

i 

i 

2 

Indicates 2 states, 1 control, 1 measurement, 

1 noise, 1 transfer function desired, no eigen- 
values, singularity of FH unknown 

21 


3 1 

Transfer function desired is for row 1, 
column 1; zeros, residues, and bode plots are 
desired; plot bounds are -1 to 1 and normal 
size plots should be done 

22 


■ 3A 

1 coefficient in the numerator of polynomial 
multiplying transfer function before plotting; 
value is 1 .0 












Table 2.3 (Continued) 


LINE 

NUMBER(S) 

GROUP 

TYPE 

CARD 

TYPE 

DESCRIPTION 

23 


3B 

1 coefficient in the denominator of polynomial 
multiplying transfer function before plotting; 
value is 1 .0 

24-25 

! 1 1 

4 

FH matrix 

26-27 


5 

GH matrix 

28 


6 

HH matrix 

29 

1 ' i 

7 

DH matrix 

30-31 


8 

TH matrix 

: 32-33 


9 

QH matrix 

34 

4 

^ 1 

Comment card indicating start of augmented 
states input 

35 


2 

No. of augmented states is 2 

36-37 , 

^ ^ 1 

3 

FA matrix 

38-39 

, : i 

^ 1 ■ 

4 

GA matrix 

40 

: 1 
i 

i 

; 5 

j 

HA matrix 

41 


6 

DA matrix 

42 


7 

VA matrix 

43 

6 

1 

Comment card indicating start of closed loop 
i nput 

44 


2 

Mo. of rows in unaugmented HC is 6, 1 transfer 
function desired, no eigenvalues, KC to be 
set to identity matrix 

45 


3 

Transfer function desired is for row 1 , 
column 1; zeros, residues and bode plots with 
bounds of -1 to 1 are desired; normal size 
plots 


56 




LINE 

NUMBER(S) 

GROUP 

TYPE 

* i 

46 

1 1 

47 

i 

■ 

48-49 


50 


51 

7 

i 

! 52 

i 

1 

j 53 


54 



m 

I 

DESCRIPTION 

3A . 

1 coefficient in numerator of polynomial 
multiplying transfer function before plotting; 
value is 1.0 

3B 

f 

1 

1 coefficient in denominator of polynomial 
multiplying transfer function before plotting; 
value is 1 .0 

1 4 

CC matrix for augmenting GC 

5 

BC matrix for augmenting GC 

1 

Comment card for additional human dynamics 
input 

1 

2 

g matrix 

I 

3 

W matrix 

4 

Vy matrix 















& 

0 

& 

o 

3 

3 



) 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

(J-l X 

^ T. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


•I'HANSfER FUNCTION PKOGRAH WITH SEVERAL OPTIONS 

(I) PLANT UVNAMICS ^COMPOTE Y/0 OR Y/M 
XDCT = fX + GU + GANMAW 

Y = HX + DU 

<2) HU.-iAN OPLRATOH DYNAMICS - COMPUTE U/Y 
SEE ROUTINE MllUMAN FOR STRUCTURE. 

(3) HUMAN OPERATOR + PLANT DYNAMICS 
COMPUTL E/R, WHERE 

E =(A*ADJL!*C) *(A2*ADJB2*C2) , FIRST ( ) FROM (I) ABOVE, 

SECOND ( ) FROM (2) ABOVE 

F =(DETD)*(DETD1) 

(4) CLOSED LOOP RESPONSE 

2DOT = P2 +GAMMAW , SEE ROUTINE MCLOSE FOR DESCRIPTION OF F AN 

Y = H2 

(5) 

C INPUT 

ALL INPUT IS PERFORMED IN INIT AND INIT2. 

TilbRE ARE ? GROUPS OP INPUT SETS. 

GROUP 1 = TITLE AND OPTIONS 
GROUP 2 = INPUT FOR PLANT DYNAMICS 

Group 3 = input for human dynamics 
GROUP 4 = input for AUGMENTED STATES 
GROUP 5 = INPUT FOR HUMAN OPRJRATOK + PLANT DYNAMICS 
GROUP 6 - INPUT FOR CLOSED LOOP RESPONSE 
GROUP ? = ADDITIONAL INPUT FOR HUMAN DYNAMICS 

THE PROGRAM LOGIC IS SUCH THAT IF THE OPTION IS SET SO THAT - 
HUMAN DYNAMICS ARE DONE, THE PROGRAM ITERATES ON GROUP 7 INPUT. 
OTHERy.ISE, IT WILL ITERATE STARTING AT GROUP 1. 

ALL MATRICES ARE READ BY ROWS ACCORDING TO FORMAT IN ARRAY IFMT. 
(SEE (IROUP 1 INPUT DESCRIPTION) 

ALL CARDS MARKED WITH * ARE CONDITIONAL INPUT. 

GROUP 1 ALWAYS INCLUDED • 

CARO fUR.'lAT DESCRIPTION 

1 (13A6) TIT'LE(13) = TITLE TO BE PRINTED AT THE TOP 

OF EACH PAGE 

2 NAMELISl- / OPTION / 

PLANT = T, TO DO PLANT DYNAMICS (DEFAULT=F) 

HUMAN = T, to do HUMAN DYNAMICS (DEFAULT=F) 

riOPLT = T, TO DO HUNAN OPERATOR + PLANT DYNAMICS 

(DEFAULT-F) 

CLOSED = T,TO DO CLOSED LOOP RESPONSE 
(DLFAULT=F) 

COVAK = T, TO DO COVARIANCES 
(DEFAULT=F) 

IFMT = 13 WORD ARRAY FOR VARIABLE INPUT FORMAT. 

ALL CARD MATRICES ARE READ KITH THIS FOR 
(DEFAULT = SEII,i.4) 

KINOUT = T, FOR MINIMAL OUTPUT 

BACKlvD = T, TO USE BACKWARD LEVERRIER IN 

COMPUTING ADJOINT EVEN IF NO. OF STATES 
IS LESS ‘IHAN 5. (AUTCMATICALLY TRIES 
•BACKlvARJD ^ ) . 

f- ■ 



n r 




c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


ON 

o 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


.. - — , - — , : \ii-j ijij iwv it&ivtii- r 

(Dlil'AllI/l-S) 

‘ DISK = IF THE PLAWT F AND G ARE TO 

BE READ OFF OF Ui-iljT 9 
XTEKC = 'r, TC ITEFATt; ON VY ANU/OW VA 
FCNTRL = T, TO S'WITCH FC MATRIX TO CONTROL ONLY 
(DEFAULT = F) 

ITERVV = NO. OF ITERATIONS FOR VY ( DEFAULT=W) 
ITEFVA = NO. OF ITERATIONS FOR VA ( DEFAULT = 0) 
RUOY = CONSTANT FOR VY (DEFAULT=.0314159) 

RHOA = CONSTANT FOR VA (DEFAULT=. B9247?) 

SAVSYS = T, IF FC, GC, HC, WBIG, VY , VA ARE TO 
BE WRITTEN TO DISK UNIT 21 
IFMT2 = 13 WORD ARRAY FOR VARIABLE FORMAT USED 
- FOR FILES USED WHEN D1SK=T OR SAVFCA=T. 

( DEFAULT=4E20.15 ) 

GROUF 2* INCLUDE ONLY IF PLANT=T 
CARD FORMAT _ DESCHIFTIOM 

1 13A6 1CMMNT(13) = COMMENT CARD DENOTING START OF 

PLANT DYNAMICS INPUT. STRICTLY 
FOR USER IDENTIFICATION. 

2 1615 NSP = NO. OF STATES 

NCP = NO. OF CONTROLS 
WMP = NO. OF MEASUREMENTS 

NOP = NO. OF NOISE 

KTFP = NO. OF SPECIFIC TRANSFER FUNCTIONS TO BE 
DONE 

JEIGPP = EIGENVALUE FLAG FOR P 

-I, COMPUTE AND SKIP ALL OTHER PLANT 
DYN. CALCULATIONS. 

0 , DO MOT COMPUTE 

1, COMPUTE AND CONTINUE 
JFSNGP “..SINGULARITY FLAG FOR F 

a, NOT KNOWN 

1, F KNOWN TO BE SINGULAR 
JENTP = ENTRY FLAG 

a, REGULAR CALCULATIONS 

1, IF TKN. FCn. was computed AND SAVED 
IN A PREVIOUS RUN, AND ONLY SPECIFIC 
TRANSFER FUNCTIONS ARE NOW DESIRED. 

3* 1615 READ NTFP SETS OF CARD TYPE 3-3B (NTFP MAY = 0) 

IRWCLP(l) = ROW MO. OP TKN. FCM . DESIRED 

(2) = COL. MO. OF TKN. FCN . DESIRED 

(3) = D, FOR ZEROES 

1, FOR ZEROES + RESIDUES 
. 2, FOR PODE PLOTS 

3, FOR EVERYTHING 

(4) = LOWER POUND OF PODE PLOT (LOG) 

(5) *= UPPER POUND OF BODE PLOT (LOG) 

(6) = fl, NORMAL SIZE BODE PLOT 

1, DOUBLE SIZE DOnP, PLOT 
3A* I5.7E10.5 INCLUDE ONLY IF IRiVCLP(3) .GE. 2 

..„NCOEPN = NO. OF COEFFICIENTS IN NUMERATOR 

OF THE N/D POLyWOMIAL WHICH MULTIPLIES 
THE TKN. FCM. REFOkE RLOTTING. 
COEFNP(I.J) = COEFFlClIiNTS OF THE NUMERATOR 

IE D ESC i: I! DING OrDI'K. 

3B* I5,?Elk).5 INCLUDE ONLY 1 F IRuCLP ( 3) .GE. 2 

MCOEFD = NO. OF COEFFICIENTS IN DENOMINATOR 
COEFDP(I.J) = DENOhINATOR COEFFICIENTS IN 
DESCt.NDlNG ORDER. 

4- 4X* IFMT F MATRIX INCLUDE ONLY IF JENTP = 0 

5- 5X* iFNT G WATrIX INCLUDE ONLY IF JEIGt’P ’.NE. '-T, JENTP‘=S?5' 



7-7X'' Ifc’.MT 
GHOUi* 3* Include 

CAhD FOi^riA'l 

1 I'JAo 

2 1615 


3* 1615 

3A* 15,?Eia.5 

3e* 15.7E1B.5 


D MATkix Include only it' jeigfp .ne. -i,jentp=u 
GAMKA Mm INCLDIi: ONLY IF jEIGPlP . Ntl. -1 , JE;nTP=0 

ONLY IF IIUMAW=T 
DESCKIPTION 

ICi'lMNl (13) = SEE CAkD 1, GHOUP 2 

NSH = NO. OF STATES 

NCH = NO. OF CONTKOLS 

NMH = NO. OF MEASOKEmENTS • 

NOH = NO. OF NOISE 

NTFH = NO. OF T1!N. FCNS. DESIRED 
JEIGFH •= EIGENVALUE FLAG. SEE CARD 2, GROUP 2 
JFSNGH = SINGULARITY FLAG FOR F. SEE CARD 2 , 
GROUP 2 

READ NTFH SETS OF CARD TYPE 3-3B (NTFH NAY = 0) 
IR'aCLH(6) = SEE CARD 3 , GROUP 2 
NCOEFN = SEE CARD 3A, GROUP 2 

COEFNIKI.J) =SEE CARD 3A, GROUP 2 
NCOEFD =SEE CARD 3B, GROUP 2 

COEFDH(I,J)=SEE CARD 3B, GROUP 2 


c 

4-4X 

IPMT 

F MATRIX 

V, 

c 

c 

c 

5-5X 

IFMT 

G MATRIX 

6— 6X 

IFHT 

il MATRIX 

c 

c 

f 

7-7X 

IFMf 

D MATRIX 

0\ 

tJ-8X 

IFHT 

GAMMA MATRIX 

1— ■ 

5-9X 

IFHT 

Q MATRIX 


Group 

4* INCLUDE 

ONLY IF HUMAN=T 

c 

CARD 

FOkMAT 

DESCRIPTION 

c 

1 

1 3A4 

ICrlMNT(13) , SEE CARD 1, GROUP 2 

c 

2 

15 

WSA = NO. OF AUGMENTED STATES 

c 

3 

iFi-IT 

FA MATRIX 

c 

4 

IFHT 

GA MATRIX 

c 

5 

IFMT 

IIA MATRIX 

c 

6 

IF.MT 

DA MATRIX 

c 

c 

c 

7* 

IFMT 

VA MATRIX, INCLUDE ONLY„IF COVAH=T 

GROUP 

5* INCLUDE 

ONLY IF hOPLT=T 

c 

CARD 

format 

DESCRIPTION 

c 

1 

1 3A6 

ICi|MN‘i (13) = SEE CARD 1, GROUP 2 

c 

c 

2 

15 

NTF3 = NO. OF TRANSFER FUNCTIONS TO 


1615 

I5,7E1D.5 

15,7l1u.5 


GROUP 6* INCLUDE 
CARD format 

1 1 3 A6 

2 1615 


READ WTF3 SETS OF CARD TYPE 3-3D (NTF3 KAY=0) 
IRwCL3(6) = SEE CARD 3 , GROUP 2 

NCOEFN = SEE CARD 3A, GROUP 2 

COEl M'KI ,J)= SEE CARD lA, GROUP 2 

NCOEfU = SEE CARD 3b, GROUP 2 

COBFD3(I,J)= SEE CARD 3D, GROUP 2 

ONLY IF CLOSED=T 
DESCRIPTION 

ICMMI'iT( 13) = SEE CARD 1, GROUP 2 
Nt-.C = NO. OF KOt.S IN HC 
NTFCiiNO. OF TKN FCNS. DESIRED 

JEIGFC = EIGENVALUE FLAG, SEE CARD 2, GROUP 2 
lUNTIlC = H = IDENTITY ‘MATRIX FLAG 

>!. u, Rir.r. PK .REA,n„TN 





lb 15 ■ 

15,7Li0.5 


j, li aIi.l hf bi;i TO ibiM ny 
NTFC SETS OF CARl; 3-l(F tFTPC i'lAY 


1Kv.CLC( 6) = SEE CAKO 3 , (iXOUE 2 
NCOEFN = SLE; CAKD 3A, GKfJUl" 2 

CCtFNC(I,J)= see: CARO 3A, CROUf 2 
NCOEft) = SEE CARD 3fi, GROUP 2 


COEEUCd ,J)^ 


SEE CARD 3fi, GROUP 
SEE CARO 3B, GROUP 


4 IfMT 
4 A IFHT 
5* lEMT 


CC MATFIX ( POH GC ) 

BC JiATRlX ( FOR GC ) 
INCLUDE ONLY IF IDNTUC = U 
H MATRIX 


GROUP ?* INCLUDE ONLY If liUMAN=T, ALL MATRICES FOR GAINS CALC. 
CARD fOk/lAT DESCRIPTICN 

I J3A6 ICHMNT(13) = SEE CARD 1, GROUP 2 


2 IFMT 


3 IFKT 


4 IFHT 


G SHALL MATRIX 


W BIG MATRIX 


VV MATRIX 


FILES USED — 

9 INPUT UNIT FOR F AND G I PLANT ) WHEN DISK=TRUE 
10 SCRATCH UNIT (AHJAYS NEEDED) 

H OUTPUT UNIT FOR PLANT DYNAMICS (TEMPORARY OR PERMANENT) 

12 OUTPUT UNIT FOR HUMAN DYNAMICS ( " " " ) 

13 OUTPUT UIUT FOR HOPLT DYNAMICS { “ " " ) 

14 OUTPUT UNIT FOR CLOSED LOOP ( " " " ) 

IMPLICIT REAL* 8 ( A-H , 0-Z ) 

LOGICAL PLANT, HUMAN, CLOSED, HOPLT 
, , MINOUT, BACKWD, DMATXP, DMATXH, DMATXC 

LOGICAL lEOF 
Logical covar 

LOGICAL SAVSYS, FCNTRL, ITERC , LITER, AITER 


C0^•f■10N 

/ 

LOGIC 

/ 

PLANT, 

HUMAN, 

HOPLT, 

CLOSED 

, HINOOT 




t 

BACKWD, 

DMATXP, DMATXH, DMATXC 

* 



9 

ITEKC, 

FCNTRL 

, AITER 



COiG'iOi'i 

/ 

SAVE 

/ 

SAVSYS 




CCMi'lOra 

/ 

uR ITS2 

/ 

lUFG , 

lUSYS 




COMhON 

/ 

MAX DIM 

/ 

NSPfiX, 

NCPMX , 

NKPMX , 

NOPHX , 

NSilHX, 1 




f 

NMHMX, 

NOHMX, 

NS 3 MX, 

NC3MX, 

NN3MX 




r 

NSCMX, 

NCCHX , 

NMCMX 



COMMON 

/ UNITS / 

lUP, lUH, 

1U3, 

lUC, JUNIT 



COMMON /VAR/ COVAR 
PLAMi DiNA.VlCS 

DIHEnSIOM FP{30,3i)) ,GP(3ii,4) ,IIP(1H,30) ,UPI1U,4) ,GAMMP(30,4) 
LI.-iEI.-SlOiV PP(31) ,IR'ACLP(r),4!)) ,G(30,8) ,NUEGNP(40) 

DIilLivSlG.I CUEF;.P(h,40) ,RDEGOP(4il) ,COEFDP(6,4ii) 

riUMAW UYi.AMiCS 

ClMEivalCN FH(15,1S) ,GH(15,4) ,HH(ll),15) ,DH(10,4) ,GAMMH(15,4) 
, uIf;Ei.:;,/UiV 0'!I(I5.J 5) ,FGAIN(15- l«) ,CGAl,N(4 ,15) ,IRWCLIIU6,16) 

— S '»'{>' " P H '.a N 0 L w .'*'{■ ( I 8 ) iitfNIi’. ^ if fcR) 'u~o DlHaav,, 

DJNLNSION CUEFDIi(6,18) ,GSML(4,<i) ,..PIG(4,4) ,VY(.lH,lti) 



AUL..iLl'.-it:D S'i'ATtS 

ri ■ 2 r ■ -> . -■■ - ■ . . ■ ■ ■■■■■ ^ 

S 5 ! . . j . . 

_ blr.L.'iSIO:>) lA(‘j,5), GA(5,4) , M(4,5) , UA(4,4), VA{4,4) 

til. \;V» Ot-i.K-"’njK + i'LAN’l IjVi.AmICs 


ol 1L.\51CN 3(bl) 16) ,COl.tN3<6^i6) 

XDhLOKie) ,CGl.l-D')(fj,16) ,Ik'.NGL3(6,16) 

CLOotC LOOf CYiMMICS 

OlflENSION FC(15, 35) ,GC(35,iy) ,HC(35,35) ,PC(36) ,XCOV(35,35) 
OIMti-SIOW NDi;Gf«C(1225) .COtTNC (6 ,1225) .^3DEGDC (1225) .COEFDC ( 6 , 1 225) 
DIMEiJSICN Iki;CLC(6,1225) DC(35,35) 

DIMEwSlON CC ( 4,4 ). HCORG ( 35 , 35) , EC (4,4 ) 

OTHER AKkAYS 

DlilEGSlON 6PLT(165) ,£IGK(35) ,EIGI(35) ,KORK(530fl) 

DIMENSION FuOUK( 35,35 ) 

DIMENSION Y0(36,122S) ,NDEGYU(1225) 

SET constants ACCOHDING TO PAKAMETEKS 

MXCO=6 

x.>:yu=36 

NSPMX=3a ' 

NCPi-iX = U 

N.-iP.'iX=l0 

NOPMX=4 

NSiliMX = 15 

NCH:-jX = 4 

HKrlNX = lU 

NOHMX=4 

NS3MX==3H 

nC3mX=S 

rJ;-i3?iX=10 

NSCMX=35 

NCCf-iX=ly 

NHCMX=35 


READ ALt, INPUT EXCEPT GSML, VY, WBIG AND 
DO SOME ITITlADIZATTON 


103 CONTINUE 

CALL I'ilT •( FP, GP, HP, DP, GAMMP, G, IRWCLP 
, NDEGNP, NDEGCP, COEFNP, COEFDP 
, Pll. GH, HH, GAMMU, OH, IKrtCLH, NDEGNH 
. , MDEGDH, COEFNh, COEFDH, DH 

, Ik‘nCL3, NDEGn3, NDEGD3, COEFN3, COEFD3 
. , HCORG, DC, FA, GA, HA, DA, VA 

, IkwCLC, NDEGMC, MDEGDC, COEFNC, COEFDC 

. , CC, bC 

. , mXCO, I EOF ) 

IF (IhOF ) GO TO 959 


PLANT DYNAMICS 

IF ( .NGl. PLANT) GO TO 230 

CALL MpLANT ( FP, GP, HP, DP, PP, YU , HXYU , IRWCLP, NDEGYO 
, NCLGNP, NOLGCP, COEFN'P, COEFDP, HXCO 
, BPLT, EIGR, EIGI, WORK, FwORK ) 


ilUMA?; L'/i.-AillCD 

■ f 



or, o ooo or, oooo ono ooo r. ooo oor, on 


^ If Co I i I. 

u ( . liUMf-.) . CO ‘lu ■’.n; 

hLAO l:. L, Abie, VY 

2b J COf/TIt-Ui: 

CALL IW1‘12 ( GSML, W3IG, V5(, JUNIT, lEOF ) 

IF ( rtOF ) GO TO 999 
26«i CONTILUL 

CUflFUTL THe CONTHOL AND FILTER GAINS ( LAy:ffDA AND K ) 

CALL MGAINS ( FH , GH , HH, OH, GSML, WBIG, VX, GAMMH 
, IVOKK, CGAIN, FGAIN ) 


COflFUTb Nfcrt F HATKICES FOR HUMAN AND CLOSED LOOP 
AND NEw H XATRIX FOR CLOSED LOOP 

CALL MFHFC ( Ft”, G, FH, GH, HH, DH, FA, GA, HA, DA, CGAIN, FGAIN 
, V(GKK, FHH, FC ,HCORG, HC ) 

CALL MHUMAN ( FHH, GIlH, HHH, DHH, HA, DA, FGAIN, CGAIN, PH 

, XU , NXYU , IRlsCLIi, NDEGYU, NDEGNH, MDEGDH 

. , COEFNH, COEFUH, MXCO, BPLT, EIGR, BIGI, WORK ) 


HU?iAN OPERATOR + PLANT DYNAMICS 
IF ( .NOT. i'OPLT) GO TO 4 00 

CALL MHOPLT ( P"l , YU , MXYU , IRWCL3, NDEGYU 

, NDEGN3, NDEGD3, COEFN3, COEPD3, MXCO 
. , EPLT, EIGR, EIGl, WORK , PP, PH ) 

CLOSED LOOP RESPONSE 

400 COtVriNUE 

IF ( .NOT. CLOSED ) GO TO 600 

CALL HCLOSE ( G, GH , GAMMP, FGAIN, FC, GC, HC, DC, CC, BC 
, PC, YU , MXYU , IRWCLC, NDEGYU 
, KDEGNC, NDEGDC, COEFNC, COEFDC, MXCO 
, EPLT, EIGR, EIGI, WORK, FlvORK ) 

IF( .NOT.COVAR) GO TO 450 

COMPUTE COVARIANCE 

CALL MCOV ( XCGV , wBIG , V Y , V A , FC ,GC , HC , WORK , FWORK ) 

CALL i-lCCST(XCOV, OH, GSML, CGAIN, RJX,RJ0,RJ, WORK) 


450 CONTINUE 

IF ( .NOT. ITERC 1 GO TO 500 

COMpU'iE nEk VY and/or VA 

CALL lilTLR { XCOV, VY, VA, ivORK, LITER, AITER ) 
IF ( LITER ) GO TO 5o0 
CiO TO 2&H 
5i*<« COtsT'lNUE 

SAVE THE SYSTEM 

IE ( .NO'l. SAVSYS ) GO TO 600 

CALL FCSAVE ( FC, GC, HC, WBIG, VY, VA ) 


t_ GE'i a CiEw •tibML, WafG" AND vx 




III. TIME VARYING SYSTEM RESPONSE PROGRAM (TVSR) 


3.1 DESCRIPTION 

The Time Varying System Response Program (TVSR) solves the 
differential equation 

i(t) = F^(t)z(t) G^(t)9(t) 

y(t) = H^(t)xCt) (3.1) 

The system matrices F , G and H are not computed in the 
program, but rather are read in from a disc file. The modified 
transfer function program described in the previous section, 
writes out these system matrices in a suitable format. Appendix 
A describes a stand-alone program which also can be used as 
generic input. 

The time varying response output can be of three types: 

(1) mean (expected) value of states/outputs 

(2) sample time histories 

(3) state/output covariance 

The system matrices are read in at the start and end point 
for some interval of time, which has been specified. These 
matrices need not have the same coefficients at the beginning 
and end. An algorithm interpolates between these matrices at 
intervals specified by the user. By this method, a linear 
piecewise solution may be obtained for the time varying system. 

An added feature of the code is that several time segments 
can be "pieced" together to make a longer time record. At the 
end of any one record, all relevant information is stored on 

preceding PAGE bL ' 
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disc. The next time segment uses this information as new start 
conditions . 

This sequence is shown in the following figure. 



^0 TIME 


Figure 3.1 

(1) Interval #1 (from tQ to tj^) requires one execution of 

the program. The system matrices (F^, must be 

previously stored on disc corresponding to t = t^ 
and t = t-, . Four interpolation intervals have been 
specified (corresponding to time points t, , t^ , 

t^ and t^ ). The state at t 2 has been specified 
for storage for initializing interval #2. (Note: t 2 

need not be equal to t^, i.e. t 2 £ t^.) 

(2) Interval #2 (from t- to t,) is now executed with 
two interpolations (t^ and t^ ) using initial 
conditions stored from the previous execution. 

The interpolation algorithm cues the value stored in the 
output vector, y, in the location INDEXV (card group 2). The 
variable, such as altitude, may be used. It is important to 
specify properly so as to ensure the variable appearing in 

the correct location in the output vector, y. 
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The structure of the code is such that the form of G 0 

c 

can be augmented from that written to disc from the transfer 
function program. The transfer function program assumes 0 is 
the form: 


0 = 


w(t) 

Vy(t) 

v,(t) 


w N(0,W) 

5 Vy N(0,Vy) 

Va NCO.V^} 



where W, and are stored on the disc. For increased 

versatility, a deterministic input of the form p(t); where 
p(t) is a polynomial function (with coefficients read into the 
TVSR program) of time. This requires augmentation of 0 of 
the form: 


p 


“ — 

r 1 G 


w 

n 



c 1 C 2 


V 



v„ 

1 G 


a 

1 c„ 



1 a 


P 

— 


* 


where now G , G and G must also be read in. 

^1 ^2 ^a 

The covariance of the augmented 0 vector is assumed to be 
of the form: 

w 0 0 0, 

0 Vy 0 0 

0 0 ° 

0 0 0 0 
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since p(t) is deterministic. This matrix is used in the state 
covariance calculation which employs OPTSYS in the solution o£ 
the steady state values. The state covariance is evaluated only 

A 

at the times designated as linear interpolation times (t^ in 
the figure) , since the system matrices only change at these 
times. The covariance matrix can be evaluated even when sample 
function option is selected, since this uses a separate algorithm 
for its calculation. 

The covariance that is printed is actually the covariance 
of y, the output. This is equivalent to the state covariance 
only when = I . 

One note of caution is required in the utilization of the 
output vector, y = H_Z. Specifically, the first NC variables 
(where NC = number of control variables - output of the 
operator) of‘ y is constrained to be u where (see Section II) 


u = -D„ A X + H„ x^ 
a a a 


Thus, y is assumed to be of the form 



1 1 



0 1 -D„ A 1 


X 

a a 

1 1 


X 

1 

a: 

o 

1 


^a 


where is read in from cards (or set t'o I internal in the 

c 

program) . 
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3.2 INPUT SPECIFICATIONS 

The program receives its input through cards and disk files. 


3.2.1 Card Input 

There are two cards required for each run. Card 1 is a 
title card for output identification. Card 2 is the option 
specification card. All succeeding input is dependent upon the 
options selected on Card 2. All matrices utilizing card input 
are read in by rows. The default format is 8E10.4 which may be 
changed by Card 2. Table 3.1 describes the required and optional 
card input. 

3.2.2 Disk Files 


The program uses four disk files: 

(1) System at start time (input only) 

(2) System at end time (input only) 

(3) Z and H*Z vectors (input/output) 

(4) Steady state and steady state output covariance 
Cinput/output) 

The user need only be concerned with the structure of 
(1) and (2). 

3. 2.2.1 File for System at Start Time 

, The unit for the system at start time is referred to as 
lUFCl and has a default value of 10. It is a sequential binary 
file. The file must contain seven logical records. All records 
must appear on the file, even if the user overrides them by 
providing card input or sets to the identity matrix. 

The format of file lUFCl is given in Table 3.2. 
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Table 3.1 
Card Input 


CARD 

^ FORMAT 

VARIABLE 

NAME 

TYPE 

DESCRIPTION 

1 

i 13A4 

ITITLE 

i 

A 

52 character title to be 
printed at the top of each 
output page 

2 

Namel i st 
/Option/ 

MEAN 

i 

L 

= T, to compute mean calcula- 
tions (Default = F) 



SAMPLE 

L 

= T, to compute the sample 
function (Default - F) 



PLOT 

L 

= T, to produce printer plots 
(Default = F) 



COVAR 

L 

1 

= T, to compute covariance 
and save on disk 
(Default = F) 



HIDENT 

L 

= T, to HI, H2 to the 

identity matrix (Default = F, 
use HI, H2 on disk files 
lUFCl, IUFC2) 

i 


WCARD 

L 

= T, to read W1 , W2 from 
cards (Default = F, use 
Wl, W2 on disk files lUFCl , 
IUFC2) 



VYCARD 

L 

; i 

= T, to read VYl, VY2 from 
cards (Default = F, use 
VYl, VY2 on disk files lUFCl , 
IUFC2) 



zcar'd 

L i 

= T, to read initial Z 
vector from cards (Default = 
F, use Z vector on disk file 
lUZ) 



VACARD 

L r 

i 

= T, to read VAl , VA2 from 
cards (Default = F, use VAl, 
VA2 on disk files lUFCl , 
IUFC2) 



lUFCl 

1 

I 

Unit No. containing FCl , 
GCl, HCl, Wl, VYl, VAl 
(Default = 10) 


A = Alphanumeric 
I =? Integer 
L - Logical 

D.P. = Double Precision 
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Table 3.1 (Continued) 


CARD 

FORMAT 

VARIABLE 
NAME , 

TYPE 

DESCRIPTION 

2 


IUFC2 

[ 

I 

Unit No. containing FC2, GC2, 
HC2, W2, VY2, VA2 (Default = 

V) 



lUZ 

I 

Unit No. containing Z, H*Z 
vectors' (Default = 12) 

'■ 


lUCOV 

I 

Unit No. containing steady 
state and steady state output 
covariance (Default = 13) 

j 


LASTZ 

I ' ■ 1 

= n, record no. of Z vector 
on file lUZ to be used as 
initial Z. (This only needs 
to be set if you do not want 
the last Z written.) 

= 0, will use last Z written 
or will start a new file if 
ZCARD = T. (Default =0) 

j 


NPRINT 

I 

Increment for printing 
(Default = 0, last Z only 
will be printed) 


j 

NSTORE 

i 

I 

Increment for saving Z and 
H*Z vectors (Default =0, 
nothing is saved) 



NINT 

I 

No. of intervals between 
TSTART and TSTOP (Default = 1) 



i NTERMS 

I 

No. of terms to use in com- 
puting PHI and PSI (transition 
matrices) (Default = 10) 



TSTART 

D.P. 

Start time (Default = 0.0) 



TSTOP 

D.P. 

Stop time (Default = 0.0) 



TDELT 

D.P. 

Delta time- between TSTART and 
TSTOP (Default =1.0) 



NCOEF 

I 

No. of coefficients in p 
matrix (j<5, Default = 3) 



ISEED 

I 

Seed for random in genera- 
tor (Default = 328765) 


A = Alphanumeric 
I = Integer 
L = Logical 

D.P. = Double Precision 
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Table 3.1 (Continued') 


CARD 




DESCRIPTION 

2 

i ’ - ' 

NPSIZE 

I 

. 

No. of rows in p matrix 
(Default=0) 



DELTAV 

D.P. 

Desired total change in 
velocity 



IMDEXV 

I ' 

Location of velocity in H*Z 
vector (Default=l) 

i 


IFMT(13) 

I 

Array containing format for 
reading matrices from cards 
(Default = 8E10.4) 

3A* 

13A4 

ICMMT 

A 

Comment card for user identi- 
fication of p matrix 

3-3n* 

IFMT 

PSML 

D.P, 

Matrix of coefficients (NPSIZE 
X NCOEF) used in computing e. 
Coefficients should be in 
decreasing order. 

*Include only if NPSIZE>0 

4A* 

13A4 

ICMMT 

A 

Comment card for user identi- 
fication of W1 matrix 

4-4n* 

IFMT 

W1 

D.P. 

Matrix (NNP x NNP) to be used 
in white noise sequence of 
sample function and/or 
covariance. 

♦Include only if WCARO = T 
and (Sample = T or 
Covar = T) . 

Hi 


nm 

A 

Same as 4A except for W2. 

fBm 

BH 

mm 

D.P. 

Same as 4-4n except for W2 


A = Alphanumeric 
I = Integer 
L = Logical 

D.P. = Double Precision 
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Table 3.1 (Continued) 


CARD 

FORMAT 

VARIABLE 

NAME 

TYPE 

DESCRIPTION 

' 6A* 

13A4 

ICMMT 

A 

Comment card for user identi- 
fication of VYl matrix 

6-6n* 

IFMT 

VYl 

D.P. ! 

Matrix (NMH x NMH) to be used 
in white noise sequence of 
sample function and/or 
covariance. 

*Include only if VYCARD = T 
and (Sample = T or Covar = T) 

Ik* 

13A4 

ICMMT 

A 

Same as 6A except for VY2 

1 7-7n* 

IFMT 

VY2 

D.P. 

Same as 6-6n except for VY2 ! 

1 8A* 

13A4 

ICMMT 

A 

Comment card for user identi- 
fication of VAl matrix 

8-8n* 

IFMT 

VA 

D.P. 

Matrix (NCA x NCA) for use in 
computing covariance. 

*Include only if VACARD = T 
and COVAR = T. 

9k* 

13A4 

ICMMT 

A 

Same as 8A except for VA2 

9-9n* 

IFMT 

VA2 

D.P. 

Same as 8-8n except for VA2 

lOA* 

13A4 

ICMMT 

A 

Comment card for user identi- 
fication of GCONEl matrix. 
*Include only if NPSIZE>0. 

10-1 On* 

IFMT 

GCONEl 

D.P. 

Matrix (NSP x NPSIZE) for 
augmenting GCl . 

llA* 

13A4 

ICMMT 

A 

Same as lOA* except for GC0NE2 

11-lln* 

IFMT 

GC0NE2 

D.P. 

Same as 10-1 On* except for 
GC2 

12A* 

13A4 

ICMMT ' 

A 

Comment card for user identi- 
fication of GCTWOl matrix. 
*Include only if MPSIZE>0. 

12-12n* 

IFMT 

GCTWOl 

D.P. 

Matrix (NSH x NPSIZE) for 
augmenti ng GCl . 

13A* 

13A4 


A 

Same as 12A* except for GCTWQ2 

13-13n* 

IFMT 

GCTW02 

j 

D.P. 

Same as 12-12n* except for GC2 


A = Alphanumeric 
I = Integer 
L = Logical 

D.P. = Double Precision 
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14A* 13A4 ICMMT A Comment card for user identi- 


fication of &CA1 matrix. 
♦Include only if NPSIZE>0. 

14- 14n* IFMT GCAl D.P. Matrix (NSA x NPSIZE) for 

I augmenting GCl 

15A* 13A4 ICMMT A Same as 14A* except for GCA2 

15- 15n* IFMT GCA2 D.P. Same as 14-14n* except for 

GC2 

:16A* 13A4 ICMMT A Comment card for user identi- 

fication of initial Z vector 

16- 16n* IFMT Z D.P. Initial Z vector (NSC long) 

♦Include only if ZCARD = T. 

;17A* ■ 13A4 ICMMT A Comment card for user identi- 

i fication of plot input 

section. ♦Include only if 

; PLOT = T. 

■17B^| I 1615 NVARZ ^ I No. of Z states to plot 

I NVARY j I No. of H^Z variables to plot 

i I i 1 

ISD I = 0, do not plot standard 

i : deviation 

[ • 

I : I ^ =1, plot standard deviation 

I : I . associated with each variable 

! ’ 1 . 

j i ISCALE I =0, use same scale for all 

- i I ^ curves 

I I =1, plot each curve on own 

i i j J. scale 

I I i KREC I = 0, plot all output 

i =1, plot only output from 

current run 

INC I = 1, plot every point 

= n, plot every nth point 

Plots are done from the out- 
put files (lUZ and lUCOV) . 

A = Alphanumeric 
I = Integer 
L = Logical 

D.P. = Double Precision 
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Table 3.1 (Concluded) 


f 


'card 

FORMAT 

VARIABLE 

NAME 

TYPE 

DESCRIPTION 

♦ — 

CO o 

. 

w—~ • r—~ 

i 

(515,5A4) 

' i 

1 : 1 
! ' i 

III 
. 1 ; ! 

i ' ! ■ 

( 

. j 

MVAR(,l-5) 

ID(i-5) 

J 

A 

j i 

; ^ 
i i 

1 i i 

Therefore, setting INC =2 
will give a plot of every 
other point on the output 
file. Maximum number of 
points that may be plotted 
is 2000.) 

List of Z and H*Z variables 
to plot 

4 character identification of 
each variable plotted 

As many sets of cards 17B*- 
17C* as desired may be 
included. Plots can only be 
done if NSTORE / 0. 


A= Alphanumeric 
I = Integer 
L = Logical 

D,P. = Double Precision 
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u 

u 

J 

u 


3. 2. 2. 2 File for System at End Time 

The unit for the system at end time is referred to as IUFC2 
and has default value of 11. Its description is identical to 
lUFCl (System at Start Time, Section 3. 2. 2.1). 


3. 2. 2. 3 File for Z and H*Z Vectors 


f The unit containing the Z and H*Z vectors is referred 

to as lUZ and has default value of 12. It is a random access 
j binary file. Its structure is defined and its records are 

written by the program. lUZ is used as an input file in two 
I cases: (1) xvhen a Z vector computed in a previous run is used 

as the initial Z vector in a later run and (2) when printer 
plots are desired. 

The structure of file lUZ is given in Table 3.3. 


Table 3.3 

Structure for File lUZ, Z and H*Z Vectors 






3.2.4 Covariance File 


The unit for the steady state and steady state output 
covariance is referred to as lUCOV and has default value of 13. 
It is a random access binary file. Its structure is defined and 
records are written by the program. It is used as an input ifile 
when printer plots with standard deviation are desired. 

The structure of file lUCOV is given in Table 3.4. 


Table 3.4 

Structure for File lUCOV, Steady State and Steady 
State Output Covariance 


^ RECORD NO, 

DESCRIPTION • 

1 

JREC = record no. of last record written on this 
file (£2000) 


NSC = no. of states (closed loop) 


NMC = no. of measurements (closed loop) 

2 -I JREC 

IREC = record no. of file lUZ 

this covariance is associated with 

' I 

YCOV = vector of the square roots of the diagonal 
elements of the steady state output covariance 
matrix (NSC long) 


ZCOV = vector of the square roots of the diagonal 
elements of the steady state covariance matrix 
(NMC long) 
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3.3 


PROGRAM FLOWCHART 


The program flowchart for the TVSR program is shown in 
Figure 3.3. 



Figure 3.3 TVSR Program Flowchart 
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3.4 TSS/360 OPERATION 


jl 

The following control cards are required to execute the 
; TVSR program on the NASA ARC TSS/360 system. 

(1) L0G0N... 

! (2) DDEF FTlOFOOl, VS, FI, DISP=OLD 

(3) DDEF FTllFOOl, VS, F2 , DISP=OLD 

(4) DDEF FT12F001, VI, Z0UT, DISP=NEW 

i (5) DDEF FT13F001, VI, C0V0UT, DISP=NEW, RET=T 

r (6) AMES IMSL 

(7) JBLB TVLIB 
C8) L0AD TVBL0CK 
(9) CALL TVflAI'N 
I (10) {data 

' (11) L0G0FF 

The description of the above execution deck, follows: 

I (1) Usual L0G0N card. 

* (2) File containing the system at start time, 

I (3) File containing the system at end time. 

i 

. i 

(4) File containing Z-vector output. 

C5) File for covariance output. 

, (6) Invoke IMSL library (for random no. generator GGN0F) . 

(7) Assign library TVLIB. 

" (8) Load in block data. 

(9) Execute program. 

(.10) Input cards. 

(11) Usual L0G0FF. 



3.5 S.VIPLE EXECUTION 


The following represents a sample execution (with descrip 
tion and output listing). 


1 3^•!FLl•- rj 

2 V OP ru g 

3 , 

4 ';sT'jP-.5, 

5 TOtlLf-. S, 

6 iA-iPLE^i. 

7 zCAPo^r, 
a .-t'd li b»J , 
g :-iCc(E7’l. 

10 *vP?.T., i‘^l , 

11 si'i'CPE^J , 

12 poorer, 

13 

14 bLLT3. , 

15 & En 0 

1 6 Pis.'iL 1*1 A r« I A 
17 ] . 

IgOCOiNSJ O.ATRTX 
1 9 1 . V’ 

2C c: . 

21 GCo.'(£2 •iATrtlx 

22 J . n 

2: d. ■> I 

!(■ GCiViOJ >1.ATRI.< 

2 £ >:■ . C 

26 b'.e 

27 GCTW02 '-lATRIX 

2 £ e.r 

29 i ) . (. 

3C cCAJ ('IATkIa 

31 0.0 

32 J.t 

32 GCA2 .-IATkIX 

34 1) . n 

39^.0 

3( IWTTIAL Z VECTOR 
3/ 0. w 
32 0 . 0 
3£ 0 . w 
4( 0.0 

41 0... 

42 0.0 

4' PLOT if:pur 
4^ 1 ] t 

49 ] 3 


' P i T . . 0 H . j h ^ 


ORIGINAL PAGE lt> 
OF POOR QUALITY 
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Table 3.5 

Description of Sample Execution 


LINE 

NUMBER(S) 

CORRESPONDING 
CARD TYPE 

DESCRIPTION 

1 

1 

Title card 

2-15 

; i'--- 1 

2 

Option card specifying start time of 0.0, 
end time of 0.5i delta time of 0.5, sample 
function, initial Z to be read in from 
cards, 1 row in the p matrix, 1 coeffi- 
cient (column) in the p matrix, every 
time point to be printed, every time point 
to be stored, printer plots to be done, 
one interval, expected velocity change is 
0.0. (Velocity change (DELTV) is meaning- 
less when no. of intervals (NINT) is 1.) 

16 

3A 

Comment card identifying p matrix 

17 

3 

p matrix 

18 

IQA 

Comment card identifying GCONEl matrix 

19-20 

10 

GCONEl matrix 

21 

llA 

Comment card identifying GC0ME2 matrix 

22-23 

11 

GC0NE2 matrix 

24 

12A 

Comment card identifying GCTWOl matrix 

25-26 

12 

GCTWOl matrix 

27 

13A 

Comment card identifying GCTW02 matrix 

28-29 

1'3 

GCTW02 matrix 

30 

14A 

Comment card identifying GCAl matrix 

31-32 

l|4 

GCAl matrix 

33 

15A 

Comment card identifying GCA2 matrix 

34-35 

15 

GCA2 matrix 

36 

16A 

Comment card identifying initial Z matrix 




Table 3.5 (Concluded) 



LINE 

NUMBER(S) 

CORRESPONDING 
CARD TYPE 

DESCRIPTION 

37-42 

16 

Initial Z vector 

43 

17A 

Comment card identifying plot input 

44 

17B 

Plot input asking for 1 Z state to be 
plotted, 1 H*Z output to be plotted, no 
standard deviation to be plotted, each curve 
to be plotted on its own scale, all output 
to be plotted, every point to be plotted 

45 

17C 

Z state 1 is to be plotted, H*Z output 3 
is to be plotted, description of Z(l) is 
IDZ, description of H*Z(3) is IDHZ 
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= iJA.Vf'tt- --- Tliit VAHKiUa PliOliilAM 
il'ItA.'J CALCULATIOij FLAG ~ F 
SA.'iPLi CALCULATIOij FLAG = T 
compute. COVAhlfifiCt FLAO = F 
I.'iFUT Wl,'-i2 FBOM CARDS FLAG = F 
IJJPUT WYl.Vya FhOfi cards flag =F ■ 
liiPUT INITIAL I FROi'l CARDS FLAG = T 
INPUT VA1.VA2 FROM CARDS FLAG = F 
PLOT FLAG = T 

SLT H TO IDSWTITI FLAG = F 
UNIT FOR FG 1. GC1 s 10 
UNIT FOR FC2, 0C2 = 11 

UNIT FOR INITIAL Z (IF DISK) AND/OR OUTPUT Z = 12 
UNIT HO. FOR COVARIANCE OUTPUT = n 
HAX. MO. RECORDS ALLOWED ON Z I/O FILE = 2000 
;-.AX. NO. RECORDS ALLOWED ON COVARIANCE FILE = 2000 


START 

TIME = 

C.OO 

STOP 

time - 

0.50 

DELTA 

TIME = 

C.G*"-. 

PRINT 

INCKEMEHT = 

1 


INCREMENT FOR SAVING Z ^ 1 

NO. OF INTERVALS FROM FC1 TO FC2 = 1 

NO. OF COEFFICIENTS IN THETA =1 

NO. OF TERMS TO COMPUTE PHI , PSI = 10 

SEED FOR RANDOM NO. GENERATOR = 328765 

FORMAT FOR READING MATRICES = (B(EIO.A) ) 

RECORD NO. OF INITIAL Z MATRIX, IF KNOWN (FOR ZCARD=F)= 0 

TOTAL VELOCITI CHANCE (DELTAV) = 0.0000 

INDEX OF H*Z (I) TO USE IN INTERPOLATION = 1 

NO. OF STATES (PLANT) = 2 

NO. OF NOISE (PLANT) = 1 

NO. OF STATES (HUMAN) = 2 

NO. OF iMEASUREfiENTS OiUMAN) = 1 

NO. OF ST A" -S (CLOSED) =6 

NO. OF INPUTS (CLOSED) =3 

NO. Or MEASUREMEHTS (CLOSED) =7 

NO. OF AUGMENTED STATES = 2 

NO. OF AUGMENTED CONTROLS = 1 

NO. OF ROUS IN PSML = 1 

SAMPLE RUN — TIME VARYING PROGRAM 


PS.ML 

1 

1 

1.00000 00 






Fcr 


1 

2 

3 

4 

5 



1 

-2.0000D 00 

0 , 0000 

0.0000 

0.0000 

0.0000 

0.0000 


2 

1 . OOOOD 00 

1 .00000-05 

1.3620D 00 

B.9443D 00 

1.80001! 01 

i.eoooo c-o 


3 

0 . 0 & 0 

9 . u ti 3 0 L 0 1 

-P.OOOOj 00 

-9.063 *.) L' 0 1 

0 . 0 V 0 0 

tj « ■. J si 


4 

C - 0 0 0 0 

1 . 3404D 01 

2 . 3f)20 D 00 

-4.51940 00 

1 • 0 .J ij *j L* 1 

1 . b 0 u u J w* 0 


5 

0 . 0 0 0 0 

0.0000 

0 . 0000 

0 . 0 0 0 0 

0 . 0 0 0 c 

1 ,0000 \j 0 0 


6 

0.0000 

0 . 0000 

-3.0710D 01 

-1.9D76D 02 

-;!.6667D Q? 

-/>,6667D 01 


oo 

o ^ 


> 

n 


*tS 


c« 


FC2 

1 

2 

3 


5 

6 

1 

-2 . OOOOD 00 

0.00.00 

0. 0000 

0.0000 

0 . 0000 

0 . 0000 

2 

1 . OOOOD 00 

1 .OOOGD-05 

3. 3820 D 00 

8, . 94 4 3D 00 _ 

1.8000D 01 

1 .8000 D_ 00 




0.0000 

0.0000 


1.00000 00 
G.OOOO 
0.0000 
0.0000 
0.0000 
0.0000 


2 

o.cooo 
0.0000 
9.06360 01 , 
1.311660 01 
0.0000 
0.0000 



O.QQOO 0.0000 0.0000 1.0000000 


3.0710D 01 -l.Sb76l) 02 -2.6b670 02 -2.66G7D 01 


3 

0.0000 
1.00000 00 
0.0000 
1.0000D 00 
0.0000 
0.0000 



1 2 3 

1 1.00000 00 0.0000 0.0000 

2 O.QOOO 0.0000 1,00000 00 

3.- 0.0000 9.06360 01 0.0000 

6 ’ 0.0000 1. 36640 01 1,00000 00 

5 0.0000 O.GOOO 0.0000 

6 0.0000 O.GOOO 0.0000 

S/iHPLE BUN Tine VARYING PROGRAM 


HC 1 

1 

2 

3 

6 

5 

6 

1 

0.0000 

0.0000 

1.36200 00 

8.96630 00 

1 .60000 01 

1.6G000 00 

2 

1.0000D GO 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

GO 3 

0.0000 

1 .00000 00 

>3.0000 

O.GOOO 

0.0000 

0.0000 

CO H 

0.0000 

0 . 0000 

1.0000D 00 

0. 0000 

0.0000 

0.0000 

6 

0.0000 

0.0000 

0.0000 

1 .00000 00 

0 . 0000 

0.0000 

6 

0.0000 

0.0000 

0.0000 

0.0000 

l.OOOOD 00 

0 , 0000 

7 

0.0000 

0.0000 

0.0000 

0 . 0000 

0.0000 

l.OOOOD GO 


ilC2 

1 

2 

3 

4 

H 

6 

1 

0.0000 

0. 0000 

1.3620D 00 

6.94430 00 

1.6000D 01 

1.60000 00 

3 

l.OOOOD 00 

0.0000 

0.0000 

0.0000 

0. 0000 

0.0000 

3 

0.0000 

1.00000 GO 

0.0000 

0.0000 

0.0000 

0.0000 

6 

0 . 0000 

0.0000 

1 .OOOCD 00 

O.DGOO 

0.0000 

0.0000 

s 

0.0000 

0.0000 

0.0000 

l.OOOOD 00 

0.0000 

0.0000 

6 

0 . 0000 

0.0000 

0.0000 

0.0000 

l.OOOOD 00 

0.0000 

7 

0.0000 

0.0000 

0.0000 

0.0000 

Ol.OOGO 

l.OOOOD 00 



1 

1 . 70001> 01 





1 1.0480D-0t 

SAKPLt fiUK TltiE VARlflNG FftOGHAM 





VA2 

1 


1 

1 . 0 * 1600-01 


GCON£1 


1 

?. 


1 

1.000CD 00 
0.0000 




0.0000 1.00000 00 




ORIGINAL PAGE IS 
OP POOR QUALITY 


5 

6 


0.0000 9.0636 

U.000& 1 

0 . 0000 0 . 

0.0000 0 . 


D 

01 

0 . 0000 

D 

01 

1 .00300 00 

0000 

0 . 0000 

0000 

0.0000 


0.0000 

0.0000 

0.0000 

0.0000 


AUGMENTED 

GC2. 







1 

2 

3 

N 



1 

1.0000D 00 

O.OGOO 

0.0000 

1.0000D 00 



2 

0 . 0000 

0.0000 

1,00000 00 

0 . 0000 



3 

0.0000 

9.0636D 01 

0.0000 

0.0000 



4 

0.0000 

1 . 3*i6JID 01 

1.000QD 00 

0.0000 



5 

0.0000 

G . 0 Q 0 0 

0.0000 

0.0000 



6 

0.0000 

0.0000 

0.0000 

0.0000 



2 0.000 

THETA 

0.00000 

0.00000 0 

.00000 0.00000 




H»2 

0,00000 

0.00000 0 

.00000 0.00000 

O.GOOOO 

0.00000 


SAMPLE RUN TIME VARXitIG PROGRAM 


INTERVAL NO. 1 

NO. OF STEPS TO BE TAKEN =10 

STARTING TIME OF THIS INTERVAL = 0,000 

End time of this interval = o.fsoo . 

INTthpOLATION CONSTANT = 0.000 

INITIAL VELOCITY = 0.000 

CUfi,TENT VELOCITY = 0.000 



1 

2 

3 

It 

5 

6 

1 

-2.0000D 00 

0.0000 

0.0000 

O.OOOO 

0.0000 

0.0000 

2 

1.00000 00 

1 .OOOOD-05 

1,3U20D 00 

8.9993D 00 

1.6000D 01 

1.8G00D 00 

3 

o.oooo 

9 .06360 01 

-2,00000 00 -9.06360 01 

0,0000 

0.0000 

N 

0.0000 

l.l'JO'tO 01 

2.36200 00 -9. ‘>1990 00 

1.80000 01 

1.60000 00 

5 

0.0000 

0.0000 

0.0000 

0,0000 

0.0000 

1.00000 00 

6 

0.0000 

0.0000 

-3.07100 01 -1.96760 02 

-2.66670 02 

-2.66670 01 



1 

2 

3 

4 



1 

1.00000 00 

0.0000 

0.0000 

1.0000D 00 



2 

0.0000 

0.0600 

1.00000 00 

0.0000 



3 

0.0000 

9.06360 01 

0,0000 

0.0000 



9 

0 . 0000 

1 .34640 01 

1.00000 00 

0.0000 



5 

0.0000 

0-0000 

0.0000 

0.0000 



6 

o.ooot. 

0.0000 

0.0000 

0 . 0000 




1 

2 

3 

4 

«> 

6 

1 

o.oooo 

0.0000 

1.36200 00 

b.9't'*3D 00 

1.60000 01 

1.60000 00 

2 

1 £ 00000 00 

0.0000 

0.0000 

0.0000 

0 . 0000 

0.0000 

3 

o.oooo 

1.00000 GO 

0.0000 

0.0000 

0.0000 

0.0000 

4 

0.0000 

O.OGOO 

1.0CC00 GO 

0.0000 

0.0000 

0.0000 

■j 

o.oooo 

0.0000 

0.0000 

1.0000D 00 

0.0000 

0.0000 

6 

0.0000 

0 . 0000 

o.oooo 

0.0030 

1.00300 00 

0.0000 

7 

0.0000 

0.0000 

0 . 0000 

0.0000 

0 . 00 00 

1.00000 00 

1 

1 

1.70000 01 







sample, run TIME VARYING PROGRAM 


0,00000 



1 

j t 




1 

J 0860D-01 

— 






1 

2 

3 

8 

5 

6 

t 

9.0488D-01 

0.0000 

0.0000 

0.0000 

oCoooo 

0.0000 

2 

5.0233D-02 

1 . 12270 00 

2.05600-02 

3.18990-03 

5.0290D-01 

6.61GOD-02 

3 

6.39030-02 

3.00520 00 

6.20930-01 

-3.0052D 00 

-2.2208D-T6 

-1.36780-17 

8 

1.7076D-02 

6.81880-01 

5. 3717D-02 

8.88030-01 

5. 02900-01 

6.6100D-02 

5 

-1 .0159C-03 

-7.5932D-02 

-2. 81210-02 

-9,66850-02 

7. 8689D-01 

2.0083D-02 

6 

-7.3901D-02 

-3.87690 00 

-9.00810-01 

-1.67620 00 

-8.89790 00 

-1 .03290-01 


PSI 



1 

2 

3 

8 

5 

6 

1 

8.75810-02 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

2 

1 .28800-03 

5. 27290-02 

9. 16160-08 

2.98880-03 

1.58 100-02 

1.87200-03 

3 

1.51210-03 

8.69270-02 

8.60690-02 

-8.69270-02 

G.OOOO 

0.0000 

8 

2.88960-08 

1,76580-02 

1.87520-03 

3.80190-02 

1.58100-02 

1.87200-03 

5 

-1.07920-05 

-1.03750-03 

-5. 18500-08 

-2. 17620-03 

8. 58330-0? 

7. 33720-08 

6 

-1 .01590-03 

-7.59320-02 

-2.81210-02 

-9.86850-02 

-2.53510-01 

00830-02 


SAHPLt nuti - — TIME VARYIHG PkOGHAM 


to 


3 0.050 THETA 

H»Z 

k 0.100 THETA 
H*2 

5 0.150 THETA 

H*Z 

6 0.200 THETA 

H»2 

7 0.250 THETA 

H»Z 

6 0. 300 THETA 

H*Z 

9 C.B*)© THETA 

H*Z 

10 O.kCO THETA 

H*Z 

11 0.J150 THETA 

H»Z 

12 0.500 THETA 

H‘Z 


0.75910E 
-.122792 
-.50'I35E 
0.1J3025E 
0 . IPg'ttiE' 
0 . 80679E 
-.20452E 
kfa61 3E 
19052E’ 
-.603‘i^E 


00- .71375E-010.66<)69E-010. 10000E 01 

&00.6 3700E-01-,27'I90E-02-.21183E 00-. 3‘<6E-0 10 .5 1 806E-020 .26297E 00 

01- .26986E-01-.3b5k5E 000.10000E 01 

0G-.n666E 00-. 13207E-01-. 12908E' 00 - . 52555E-0 1 0 . 22 9 35E-0 1 0 . 37« 38E 00 
010.60853E-010.38505E 000.10000E 01 

000. 36300E-020.9A538E-010. 185832 000 . 58606E-0 1 0 . 2869 2E-0 1 - . 27 M90E 00 
01-.5k3bOE-010.96275E 000. 10000E' 01 

00- . 96kk8E-010.6H299E-01-,51951E-010.5l237E-010.it07A2E-02-.52553E 00 

01- .'('I962E-01-. '(56'I5E OOO.IOOOOE 01 

00-.85097E-010.25396E-01-. 13176E CO- . 658 1 5E-02- . 1 23 1 8E-0 1 - . 78255E-0 1 


0.99293E 01-.Gk3C8E-020.90ii3bE OOO.IOOOOE 01 

0 . lG0k7L 00 0 . 2051 3E OOO . vj 26 3 E-0 1 - . 29*i6 92-0 1 0 . 1 A 337E-0 1 - . 9 1299E-020 . 98577 E-0 1 
0.15kd9E 01-.99591E-020.24605E' OOO.IOOOOE 01 

0.11299E 000.3O6O9E GOO . 7 1 H09E-0 10 . 8 1 770E-0 1 0 . 88621 E-G 1 - . o6702E-02- . 108 32E 00 
0.39868E 010.86170E-010.58610E OOO.IOOOOE 01 

-. 816682-020. 513C6E 000.12757E 000.26676E 000 . 1 2608E. 00- . 25 1 9 1 E-0 1 - . 60662E GO 
-.688932 010.22516E-010.31837E OOO.IOOOOE 01 

-.98363E 000.20305E 000.13615E 000.33639E 000.18671E 00-. 68 S67E-0 1 - . 8658 8 E 00 
0.17870E 01-. 667222-010, 86781E OOO.IOOOOE 01 

-.151382 010.316382 000.10207E 000 . 65 1 6 1 E-O 1 0 , 7 11 662-0 1 - . 905 3 1 E-O 1 - . 279 36E 00 


PLOT HO. 1 

fiO, OK VARIAHLES TO BE PLOTTED (Z) = 1 
•MO. Or VARIABLES TO BE PLOTTED (ii»Z) = 1 
STANDARD DEVIATIOfl PLOT FLAG = 0 

SCALE KLAG = 1 

ENTIRE FILE PLOT FLAG = 0 

PLOT Ii.CitEMENT = 1 

Z STA'’’ES To Be. HLOTTEC = 1 

il*Z STATE.i ■■'O BE plotted = 3 

0.000 0.000000 0.000000 
0.0500. 8370 08i-01-.278900E— 02 
0. ICO-. 1 16659L GO- . 1 32070E-01 
0.1500. 3 d 3008E-020. 8853792-01 
0.200-. 8 68 48 IE-0 10.68299 12-01 
0.250-. 8509 7 32-010.25395 7 E-O 1 
0. 3000. 205 1262 000 . 3926 30E-0 1 

f' 



r~j 


o o 





o o o 




. ‘lOio .‘j 1 '-(ObOt 000 . i;'V‘)bOX uU 

, 4 ‘i 00 . 2030 '> 1 t 000 ,n 6 ' 50 E 00 

.500G.316^'i6£ 000. 1G2b71t 00 

STATE. -'H.OOt-OI -2.C0E-01 -l.OOE-01 0.00 1.00E-0.1 ?.G0t-01 H.OOE-01 <l.0&c-01 O.OOE-01 6.0Gc-01 7.00E-01 

liiZ + - + + - — — 4 . -.4 — + 

STATE -;>.00c.-0;» 0.00 2.00E-0;> f(.00t’-GR 6.00t-0^ 3.001-0? 1.001-01 1.?01-01 1.401-01 1.601-01 

ICHZ + + + 4 - 

Tlrit 




0 . 0 G 0 G 0 



0 


A 









0. 05000 


_ 

B' 



A 








G. 10000 


ii 

A ' 











0. 15000 





A 

B 








0.20000 




A 

' 


B 







0.25000 



'A 


' B 









0.30000 





' 

B 

A 







0.35000 







B 

A 






0.40000 





* 





A B 



- 

0.45000 





* 


A 




B 



0.50000 





* 



A 

B 








<~o 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


fhOGKAH TO SOLVE 

ZDOT = t'C*Z + GC*THETA 

THE SOLUTION IS GIVEN BY 

Z(K+1) = BHI^ZiK) + PS1*GC*THETA 
mHLHE nil = STATE TkANSITTON MATRIX 
PSI = INPUT TKANSITTON MATRIX 


vl 

THETA = a , FOR MEAN CALCULATIONS 

PSML 

W , FOR SAMPLE FUNCTION WHITE NOISE SEQUENCE 

THETA = VV GENERATED FROM N ( 0 , W) , N ( 0 , VY ) & N(0,VA) 

VA 

PSML 


INPUT — ( * CARDS ARE CONDITIONAL ) 

Card mo. format description 

1 13A6 TITLE = TITLE TO BE PRINTED AT TOP OF PAGE 

2 NAMELIST /OPTION/ 

MEAN = T, FOR MEAN CALCULATIONS 
= F, DO NOT COMPUTE MEAN 
( DEFAULT-F ) 

SAMPLE = T, TO COMPOTE THE SAMPLE FUNCTION 
= F, DO NOT COMPUTE 
( DEFAULT = F ) 

PLOT = T, TO PLOT ON PRINTER 
= F, DO NOT PLOT 
( DEFAULT = F ) 

WCARD = T , TO READ IN V<1,W2 FROM CARDS 

= F, USE Wl,H2 MATRICES STORED ON DISK 
( DEFAULT = F) 

VYCARD = T, TO READ VY1,VY2 FROM CARDS 

= F, TO USE VYl, VY2 STORED ON DISK 
( DEFAULT=F ) 

ZCARD = T, IP INITIAL Z VECTOR IS CARD INPUT 
= F, IF IMTIAL Z VECTOR IS OM DISK 
( DEFAULT=F ) 

VACARD = T, TO READ VA1.VA2 FROM CARDS 
= P, USE VA1,VA2 STORED ON DISK 
( DEFAULT = F ) 

NPSIZE = NO. OF KDa'S JN PSML ( DEFAULTER ) 

COVAR = T, TO COliPUTE COVARIANCE AND SAVE 
= F, DO MOT COMPUTE COVARIANCE 
(DEFAULT=F ) 

lUFCl = UNIT NO. OF FCl, GCl , Wl , VYl, PI, HCl 
( DEFAULT=10 ) 

IUFC2 = UNIT WO. FOR FC2, GC2, W2, VY2, P2, HC2 
( DEFAULT=IJ ) 

lUZ = unit wo. FOR INPUT/OUTPUT Z VECTORS 
( DEFAULT=12 ) 

lUCOV = UNIT WO. FOR COVARIANCE OUTPUT, IF ANY 
^ t ' « ' • r -'AUL' '" ) ■ ' ■ ' 

‘ ‘ ” * L/{STZ’= record no. of initial Z 'fo BE USED ON 

Fil,r. 'IH)S OMi.Y Mrpps 'in PK <;i'T if 







n n n n d n n n r 


NPRINT 


NSTOKE 


NTERMS = 


TSTART = 


•rSTOP = 


TDELT = 


NCOEf = 


ISEEO = 


HIDENT = 


DtXTAV 

IMDEXV 


= U, ImIEL use LASl Z rtklTTEN Ok wIEL STAUT 
A E1L|' fk ZCAS"''-'l'WUl:. : _ ; 

( ' bui AUL'i f . 

= IWCKEHENT FOK PRINTING 

( 0EE'AULT=0, ONLY LAST Z PHIMTED) 

= IMCKEHENT FOR STORING 2 VECTOR, 

EVERY NSTORE VECTOR WILL BE SAVED 
( CEFAULT=a, NOTHING IS SAVED) 

NO. OF INTERVALS BETWEEN TSTART, TSTOP 
(FOR INTERPOLATION EETv\'EEN FC 'S ET. AL.) 

( DEFAULT=1 ) 

= NO. OF TERMS TO USE IM COMPUTING 
PHI, PSl 
( DEFAULT = 10 ) 

= START TIME 

( DEFAULT = H.0) 

= END TIME 

( DEFAULT=0.0) 

= DELTA 'iTNE BETWEEN TSTART, TSTOP 
( DEFAULT=1.0) 

= NO. OF COEFFICIENTS IN PSML MATRIX 
( COL. DIMENSION OF PSML ) 

( DEFAULT=3 ) 

= SEED FOR RANDOM NO. GENERATOR 
( DEFAULT=326?65 ) 

= FORMAT FOK CARD MATRIX INPUT 
( DEFAULT = ecia.5 ) 

= T, SET H TO IDENTITY MATRIX 
= F, USE ORIGIINAL II 
( DEFAULT=F ) 

= DESIRED TOTAL CHANGE IN VELOCITY 
= LOCATION OF VELOCITY IN OUTPUT VECTOR 
H*Z {DEFAULT=1) 


3* 13A6 ICMMT = COMMENT CARD TO IDENTIFY PSML MATRIX 

FOR user's PURPOSE ONLY 
INCLUDE ONLY IF NPSIZE .GT. 0 

4-4N* IFMT PSML MATRIX, READ IN BY ROWS 

NO. OF ROWS OF PSML = NPSIZE 
INCLUDE ONLY IF NPSIZE .GT. I) 


5A* 13A6 


5-5N* IFKT 


ICMMT = COMMENT CARD FOR USER IDENTIFICATION 
OF Wl MATRIX INPUT. 

INCLUDE ONLY IF ISAHPLE=T OR COVAk=T) 
INCLUDE ONLY IF MEAN f FALSE OK COVAK=TRUE 
AND WCARD=T 

Wl MATRIX 


6A* 13A6 ICHMT = SAME AS 5A* EXCEPT FOR W2 

6-6N* IFHT SAME AS 5* EXCEPT FOR W2 


7A* 13A6 

7-7N* IFMT 


ICMMT = SAME AS 5A* EXCEPT FOR VYi 
SANE AS 5* EXCEPT FOR VYI 


faA» 13A6 
E-BN* IFMT 
9A* 13A6 
9-9N* IFMT 


ICMMT = SAME AS 5A* EXCEPT FOR VY2 
SAME AS 5* EXCEPT FOR VY2 
SAME A 5A* EXCEPT FOR VAl 
SAME AS 5* EXCEPT FOR VAl 


IHA* 13A6 SAME AS 5A* EXCEPT FOR VA2 

lO-lBN* SAME AS 5* EXCEPT FOR VA2 


llA* 13A6 
11-llN* IFMT 


COMMENT CARD FOR START OF GCONEl 
GCONEl MATRIX 



c 

12A* 

13A6 

COMMENT CARD FOR 

START 

OF 

GCONE2 

c 

c 

c 

12-12N* 

IFMT 

GG0NE2 MATRIX 




13A» 

13A6 

COMMENT CARD FOR 

START 

OF 

CCTWOl 

c 

13-13N* 

IFMT 

GCTWOl MATRIX 




c 

14 A* 

13A6 

COMMENT CARD FOR 

START 

OF 

GCTW02 

c 

14-14N* 

IFMT 

GCTW02 MATRIX 




C 

c 

15A* 

13A6 

COMMENT CARD FOR 

START 

OF 

GCAl 

c 

15-15N* 

IFMT 

GCAl MATRIX 




c 

16A* 

13A6 

COMMENT CARD FOR 

START 

OF 

GCA2 

c 

16-16W* 

IFMT 

GCA2 MATRIX 





c 

C X7A* 13A6 

C 

C 17-17W* 

C 
C 

C 18A* 13A6 

C 
C 

C 18B* 1615 

C 

c 
c 
c 
c 
c 
c 
c 
c 

'-O C 

c 

C 18C (515, 5A4) 

C 
C 
C 
C 

C GENERAL NOTES — 

C (1) ALL MATRICES ARE READ IN BY ROWS FOR CARD INPUT 
C 

c (2) CARDS :j-S7n are included only if meam=t or SAMFLE=T 


C 

3-' 

THRU 

4N* 

ARE 

USED 

ONLY 

It NPSIZE .GT. 0 

C 

5A* 

THRU 

6N* 

ARE 

USED 

ONLY 

IF sample=true ok covar=true 

C 







AMD WCARD=TRUE 

C 

7A* 

THRU 

8N* 

ARE 

USED 

ONLY 

IF SAriPLE=TRUE OK COVAR=TRUE 

C 







AND VYCARD=TRUE 

C 

9A* 

THRU 

10N* 

ARE 

USED 

ONLY 

IF COVAR=TRUE, VACARD=TRUE, 

c 







AMD NPSIZE .GT. 0 

c 

llA* 

THRU 

16N* 

ARE 

USED 

ONLY 

IF NPSIZE .GT. 0 

c 

17A* 

THRU 

17 W* 

ARE 

USED 

ORLY 

IF ZCARO = TRUE 

c 

18A* 

THRU 

IHC* 

ARE 

Ui.r:D 

ONLY 

IF PLOT = TRUE 

c 

c 

(3) UNIT USED- 







c 

5 - CARD READ 






c 

6 - PRINT 







c 

IH - INPUT 

FILE 

CONTAINING FCl 

, GCl 

, HCl , ivl, VYl, VAl 

c 

11 - 

' ' 


. ' ' 

FC2 

:, GC2 

, HC2, W2, VY2, VA2 

c 

12 - IN/OUT 

' ' ' 



Z VECTORS 

c 

13 - OUTPUT 

FILE 

FOR 

COVARIANCE, RANDOM ACCESS 

c 

.AKAR^^_3 NE- 

■ . '1 = 

4,.... 

^4,N. 

. . ,NM*' . , 

NCO ... GCA- 


C PARAMETER UPL'i-lriO 


ICHMT = COMMENT CARD FOR USER IDENTIFICATION 
OF INITIAL Z VECTOR. INCLUDE ONLY IF ZCARD=TRUE. 
AS MANY CARDS AS NEEDED FOR Z VECTOR INITIAL 
VALUE. INCLUDE ONLY IF ZCARD=TRUE. 

COMMENT CARD FOR IDENTIFICATION OF PLOT INPUT 
TOTAL NO. OF VARIABLES THAT MAY BE SPECIFIED 
MUST BE .LE. 5 

NVARZ = NO. OF Z STATES TO PLOT 
NVARY = NO. OF 1I*Z VARIABLES TO PLOT 

(DO NOT INCLUDE STD. DEV. CURVES IN NO.) 
ISD = 0, DO NOT PLOT STD. DEV. 

= 1, PLOT STD. DEV. ASSOCIATED WITH EACH 
VARIABLE. 

ISCALE = 0, USE SAME SCALE FOR ALL CURVES 

= 1, USE DIFFERENT SCALE FOR EACH CURVE ' 
KREC = 0, PLOT ALL OUTPUT 

= 1, PLOT ONLY THE CURRENT RON 
INC = 1, PLOT EVERY POINT 

= N, PLOT EVERY NTH POINT 
HVARZ(5) = LIST OF Z OR H*Z STATES TO PLOT 
( COLS. 1-25 ) 

Y1D(5J = 4 CHARACTER DESCRIPTION OF EACH 
VARIABLE PLOTTED (COLS. 26-45) 






.\ZKI.C=it/y , NCUiU-,C=J i;.) 


1 (; ITS "'Hfc* A i. A”il .•'■a — p- ) * ’ " 

ht'AL fLV 

LOGICAL MEAN, wCAKU, VVCAKt), 2 CAKD, COVA«, VACAKD 
, SAMfLE, ifLOT, UlDEI-ri* 


COMMON / MA.XDIM / MXNS, MXNl , MXHM, MXNNP, MXK'MJ*, MXNCO, KXNGCA 
. , MXPLT, MXZ, MXCOV 

COMMON / LOGICS / MEAN, COVAH, wCAKD, VYCARD, ZCARD, VACARD 
, SAMPLE. PLOT, HI DENT 
COMMON / UNITS / lUfCl, IUEC2, lUZ, lUCOV 


DIMENSION 


FC1{ NS. NS ), I’C2( NS, MS ), fCC ( NS, MS ) 

GC1( NS,NI ), GC2( NS, MI ), GCC ( NS,NI ) 

•,•.*1 (NKP.NNP), W2 (NNP.NNP), WC (NNP.NNP) 

VYl (NMH.NHH) , V Y2 ( NMH , NMH ) , VVC(NMH,NMn) 

HC1( NM.H3 ), HC2( NM.MS ) , HCC ( HM,NS ) 

PI( NGCA. NGCA ), P2( MGCA , NGCA ), PC ( NGCA , NGCA ) 
PHI( NS, NS ), PSK NS, MS ), THETA( NI,2) 

2( NS,2) , PSML( NGCA, MCO ) 

ZCOV( MS, NS ), YCOV( NM.NM ), PS1GC( NS.NI J 
PLT<NPLT.ll) 

V\OKK ( ivRK ) 


DIMENSION tCi( 35, 3S ), FC2 ( 35,35 ), tCC ( 35,35 ) 

, GC1( 35,19), GC2( 35,19 ), GCC ( 35,19 ) - 
, W1 ( 4 , 4 ) , W2 ( 4 , 4 ) , WX) < 4 , 4 ) - 
, VYl ), VV2(lll,lb ), VYC( iw.ia ) - 


, HC1 (35,35 ), HC2 (35,35 ), liCCI 35,35 ) - 

, PH1 ( 35,35 ), PSI (35,35), T1IETA( 19,2 ) - 
, Z( 35,2 ), PSHL ( 5,5 ) - 

, ZCOV ( 35,35 ), VCOV ( 35,35 ), PSIGC ( 35,19 ) - 

, PLT ( 2;)DU, 11 ) 

, NORK ( ) 

DIMENSION fCCO\/( 35,35 ), GCORGl ( 3 5 . 19 ) , GCORG2 ( 35 , 19) 

, GCONEl(3t),5) , GCOME2 (30,5) 

, GCTk\01(30,5) , GCnc02(3u,5) , GCAl(5,5), GCA2(5,5) 

, VA1(4,4), VA2(4,'4), VAC (4, 4) 


EXTERNAL TIIETAM, THETAS 
CONSTANTS BASED ON PARAMETERS 


C 

MXNS = 

MS 

c 

HXNI = 

Ml 

c 

MXNM - 

NM 

c 

MXNNP 

= NNP 

c 

nXNMH 

= NMH 

c 

HXNCO= 

NCO 

c 

HXNGCA 

= NGCA 

c 

nxz = 

HZ EEC 

c 

MXCOV 

= NCO R EC 

c 

MX PLT 

= NPLT 

c 



c 

MXNS = 

35 


MXNl = 

19 


I'iXt-iM = 

35 


MXNNP 

= 4 


MXt.Mrt 

= 10 


MXIvCO 

= 5 


MX NGCA 

= 5 


MXZ = 2000 
MXCOV = 2000 


m 



i I A V Lj 


BEAD THE OPTIONS, SET ALL DIHEilSlONS, CHECK SOME VALUES 

lai) CONTINUE 
CALL IlillL 

INPUT ALL MATPICES 

IF ( .NOT. MEAN .AND. .NOT. SAMPLE ) GO TO 300 
CALL INPUT ( PSML, FCl, FC2, GCOKGl, GCCHG2, riCI , IJC2 
, Wl, w2, VYl, VX2, 2, VAX, VA2 
. GCONEI, GCONE2, GCiV^Ol, GCTW02, GCAl, GCA2 ) 

FOKM THE AUGMENTED GC AND HC MATRICES 

CALL FOKMGC { GCORGl , GCOWEl , GC'lV.Ol, GCAl, GCl ) 

CALL FORIIGC ( GCORG2, GCONE2, GCTW02, GCA2, GC2 ) 

IF ( .NOT. MEAN ) GO TO 200 

COMPUTE Z FOR MEAN CALCULATIONS 

CALL MINTEG ( FCl, FC2, FCC, GCl, GC2, GCC, HCl, HC2, HCC 
. Wl, W2, KC. VYl, VY2, VYC, VAX, VA2, VAC 

• . Z, PSML, PHI, PSI, PSIGC, THETA, ZCOV, YCOV, WORK 
, FCCOV, THETAM ) 

GO TO 3b0 

COMPUTE 2 FOR SAMPLE FUNCTION 
20U CONTINUE 

CALL KINTEG ( FCl, FC2, FCC, GCl, GC2, GCC, HCl, HC2, HCC 

• . Wl, «2, WC, VYl, VY2, VYC, VAl, VA2, VAC 

. , 2, PSML, PHI, PSI, PSIGC, THETA, ZCOV, YCOV, WORK 

. , FCCOV, THETAS ) 

PLOTTING 

3U0 CONTINUE 

IF ( .NOT. PLOT ) GO TO 9 SR 

CALL MPLOT ( PLT, 2, WORK, ZCOV,_YCOV ) 


999 CONTINUE 
STOP 



IV. OPTIMAL SIMULATION PROGRAM (TVOPT) 

I The time varying optimal simulation program (TVOPT) , as 
previously discussed, utilizes an exact closed loop simulation. 
That is, the time delay T, in the perceived display and Kalman 

filter predictor equations, is explicitly accounted for. (The 

1 

program TVSR relies on a lumped linear approximation to this 
delay time.) The program TVOPT is more restrictive than TVOPT 
because it does not compute covariance terms but is limited to 
sample time histories. 

4.1 DESCRIPTION 

The system to be simulated is given by: 
plant 

x(t) = Fx(t) + Gu(t) + Tw(t) 
y(t) = Hx(t) + Vy(t) 

Vp(t) = y(t-T) 

estimator 

x(t-Tjt-T) = Fx(t-Tlt-T) + Gu(t-T) 

+ K[y^(t) - Hx(t-T|t-T)] 

predictor 

^(tjt-T) = Fx(tjt-T) + Gu(t) 

+ e^'^K[ypH) - Hx(t-Tlt-T)] 

control 

u(t) = -Xx(t I t-T) 
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Defining y(t) = xCt|t-T) and p(t+T) = q(t) = x(t|t), then 
z(t) = z(t) + 0(t) 


where 


H 


n 


z'f(t) = [x'^(t). Y^'ct), q^'ct)] 


»: i 

li 


F = 


F -GX 0 

0 * F-GX 0 

KH -GX F-KH 


1 -! 


G = 


G r 0 0 

GOO e^’^K 


G 0 K 0 


e^ct) = W^t), Vy^ct), bVt-T)] 

b(t) = V^(t) + H[xH) - q(t)] 

The displayed variable (input to the operator) is given by 


XpCt-T) where 


yp(t) = [H 0 0]z(t) + Vy(t) 


and the operator’s output, given by 


= [0 -X 0]z(t) + v^Ct) 


' I 

m» 


A 
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It is these last two equations which are computed as 
written to an output disc file (/pCt-T) and u^^Ct)). 

The solution for z(t) is obtained from 
z(k+l) = (^> 2 ^ z(k3 + 2 ® Ck) 

i\rhere 

= exp A ; A = integration step size 

4>2 = ]. <J>i dt 
0 

The time delay, T, is restricted to be integer values of 
the integration step size, i.e. T = NA; N is input. The only 
user option is mean values of the states. This is implemented 
in the same way as the sample function case but with the a - 
values on all noise sources sbt to zero. The user must specify 
the matrices F, G, H, T, Q (in the optimal cost), W (= cov w) , 
Vy. (= cov Vy.) , Vjjj (= cov Vjjj) , g (in the optimal cost), and 
X (t = 0) . 

4.2 INPUT SPECIFICATIONS ' 

All input is accomplished through cards. All matrices are 
read in by rows in format (qElO.4). Table 4.1 describes required 
and optional card input. 
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Table 4.1 

Card Input Specifications 



(13A4) ITITLE(13) = title to be printed at the top of each 
output page 


2 NAMELIST MEAN=T, to compute mean calculations (Default=F) 

/TVOPT/ 

SAMPLE=T, to compute sample function (Default=F) 

STAT=T, to compute mean and standard deviation for 
W, Vy, u, y and Z(l-NS). (Default=F) 

PLJ3T=T, to get printer plots of z, u, and/or y 
(Default=F) 

DELTA = delta time (Default=l .0) 

N = time delay (Default=l) 

MPTS = total no. of points to compute (Default=10) 

I SAVE = increment for saving points on units 10 and 

II (Default=l) 

NS = no. of states (no default) 

1 

NC = no. of controls (no default) 

N0 = no. of outputs (no default) 

NM = no. of measurements (no default) 

NTERMS = no. of terms in expansion series in 
computing (p and ip (routine DISC) (Default=10) 

ISEED = seed for random no. generator (Default= 
1487621) 

3- 3N (8E10.4) F matrix (NS x NS) 

4- 4N (8E10.4) G matrix (NS x MC) 

5- 5N (8E10.4) H matrix (NM X NS) 

6- 6N (8E10.4) r matrix (NS x N0) 
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Table 4.1 (Continued) 


PAGE IS 
i- QUALITY 


; CARD NO. 

1 

FORMAT 

DESCRIPTION 

7-7N 

(8E10.4) 

Q matrix (NS x NS) 

8-8N 

(8E10.4) 

W matrix (N0 x N0) 

9-9N 

(8E10.4) 

V matrix (NM x NM) 

1 10-1 ON 

(8E10.4) 

matrix (NC x NC) 

11-llN 

(8E10.4) 

g matrix (NC x NC) 

12-12N 

(8E10.4) 

XO matrix (NS x 1 ) 

13* 

m 

NAMELIST 

/PLOTS/ 

NZPL0T = no. of z variables to plot (max. of 5) 
(Default=0) 

NUPL0T = no. of u variables to plot (max. of 5) 
(Default=0) 

NYPL0T = no. of y variables to plot (max. of 5) 
(Default=0) 

NZ(l-5) = list of z variables to plot 

NU(l-5) = list of u variables to plot 

NY (1-5) = list of y variables to plot 

INC = plot increment (this will plot every n— 
point stored on unit 10 and/or unit 11. I.e., if 

ISAVE=2 and INC=2, every 4HL point computed will be 
plotted.) (Default=l) 

INDSCZ=T, to independently scale each z variable 
(Default=F) 

INDSCU=T, to independently scale each u variable 
(Default=F) 

INDSCY=T, to independently scale each y variable 
(Default=F) 

PRINT=T, to echo the data to be plotted (Default=F) 
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CARD NO 


FORMAT 


DESCRIPTION 


14-14N 4A4 YLAB(4, 1-NZPLOT+NVPLOT+NYPLOT) = label of up to 16 f! | 

characters. One card per variable in order of z ? 

variables, u variables and y variables. j 

ii» 'll I I III I I II I ■■ I ■ ■ I I ■ ■■ ■ I II II ~ ■ I I i 

^ t ? 

*Cards 13*— 14* are only required when PL0T=T. As many sets of cards 13- 
14 as desired may be stacked. 
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4.3 OUTPUT FILES 


The program uses 2 output files. Unit 10 contains z-vector 
output and unit 11 contains u-vector and y-vector output. Tables 
4.2 and 4.3 describe these files. 


Table 4.2 

Z-Vector Output (Unit 10) 


RECORD NO. 

DESCRIPTION** 

1 

TIME, z(l), 2(2), ..., z(MSC)*** 

2 

• 

• 

• 

• 

• 

• 

MSAVED* 

• 

• 

• 


*NSAVED = NPTS/ISAVE 

**A11 output variables are double precision 
***NSC = 3** no. of states 


Table 4.3 

U-Vector and Y-Vector Output (Unit 11) 


RECORD NO. 

DESCRIPTION** 

1 

TIME, u(l), u(2) u(NC),*** y(l), y(2), . . . , y(NM)**** 

2 

• 

• 

• 

NSAVED* 

• 

• 

• 

• 


*NSAVED = NPTS/ISAVE 

**A11 output variables are double precision 

***NC = no. of controls 

= no. of measurements 
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4.4 PROGRAM FLOWCHART OF TVOPT 

The program flowchart for the TVOPT program is shown in 
Figure 4.1. 



f I 
11 






INPUT OPTIONS 
ANO MATRICES 


COMPUTE ANO <>2 


OHIGINAL page is 

<-*E POOR QUALITY 



COMPUTE CONTROL AND 
FILTER GAINS 



I" 

COMPUTE 


! 

* 

r" 

FORM FC ANO GC 


i 


;^AN-TRUE?; 

[NO 


COMPUTE z. u, y FOR 
MEAN CALCULATIONS 


WLE»TRU^ 

INO 


COMPUTE z, u. y FOR 
SAMPLE FUNCTION 


COMPUTE MEANS ANO 
STAnOARD DEVIATIONS 


^TAT-TRUEpi#. 

|NO 

'^LOT*TRUE?^ 

^ES 


REAO PLOT INPUT 


PRODUCE PRINTER PLOTS 
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4.5 TSS/360 OPERATION 

The following control cards are required to execute the 
TVOPT program on the NASA ARC TSS/360 system. 


(1) 

L0G0N. . . 


(2) 

DDEF FTlOFOOl, 

VS, Z0UT, DISP=NEW, RET=T 

(3) 

DDEF FTllFOOl, 

VS, UY0UT, D ISP* NEW 

(4) 

AMES IMSL 


(5) 

JBLB TV0PT.LIB 


(6) 

CALL MAIN$$ 


(7) 

{input data 


(8) 

L0G0FF 


The 

description of 

the above execution deck follows 


(1) Usual L0G0N card. 

(2) Data definition card for unit 10, z-vector output. 

The DISP=NEW parameter indicates the file named Z0UT 
does not exist. The RET=T parameter CRETAIN=TEMPORARY) 
means the file will automatically be erased at the end 
of the job. 

(3) Data definition card for unit 11, u-vector and y- 
vector output. Since the RET=T parameter does not 
appear on this card, the file UY0UT will be permanent. 
The user will have to use the ERASE command in a later 
run to delete it. 

(4) Invokes the IMSL library. Routine GGN0F (random no. 
generator) is used from IMSL for the sample function. 

(5) Specifies that the library TV0PT.LIB should be searched 
for all routines first. 

(6) Execute the program. 

(7) Input data. 

C8) Usual L0G0FF card. 


108 


ORIGINAL page IS 

OF POOH QUAT/tlY 


4.6 SAMPLE EXECUTION 


The following represents a sample execution (with descrip' 
tion and output listing) . 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

n 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 
23 


51 

52i 

53: 

54 


'If’ TD \.;C Tv 

CL03£0*T, 

SAVaYS •T , 

V. oVAK*"r . 
iS.iO 

PC?, PLAOT CilNA.'ITCS 
2 1 ] 1 
- 2.0 0 . i> 

1.0 l.i.£.'-0S 

i'. 

1 . 1 ! 

I) . 0 
d.«) 
l.i) 

>J. c 

IWPUT POR aO.IAN CWAMJCS 


1.0 


2 

1 

] 

1 


1 
1 

1.2 

1.0 


24 

-2 , tl 


.0 


25 

i.o 


1.0E 

-0 5 

26 

0.0 




27 

1.3 




28 

0.0 ! 


• III 


29 

0.0 




30 

1.0 




31 

0.0 




32 

O.G 

■3 

t V 


33 

0.3 

] 

• d 


34 

AUGMENTS 

D STATES 

INPUT 

35 

2 




36 

e.e 

1 

. 0 


37 

-266 . 67 

- 

26 . 667 


38 

0.0 




39 

44.444 




40 

13. 

} 

.3 


41 

-2.0 




42 

0.104S 




43 

CLOSED LOOP 

INPUT 


44 

6 

1 

0 

1 

45 

1 

1 

3 

-1 

46 

1 J 

.0 



47 

1 1 

. 3 



48 

0.0 




49 

0.0 




50 

1.0 





ADOITTOWAiL. HUMAN INPUT 

.05 

17, 

. P3123 
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V. IDENTIFICATION PROGRAM (SCIDNT) 


The identification algorithm delivered to NASA ARC is a 
modified version of an aircraft identification code commonly 
referred to as SCIDNT. The basic modification utilizes the 
output error method to estimate the parameters associated with 
the human operator model. Specifically, the parameters can 
include any combinations of the parameters associated with the 
optimal cost function as well as the parameters associated with 
the Kalman filter/predictor (e.g. the process noise/measurement 
noise covariances as well as the prediction time) . 


5.1 DESCRIPTION 

The identification algorithm assumes a structure for the 
human operator model of the form: 

x(t|t-T) = F* x(tlt-T) + G* y^Ct-T) 


u(t) = -X x(t t-T) + Vjjj(t) 


(5.1) 


where 


F* = (I + e^"^ KHT) (F - GX - e^"^ KH) 


G* * (I + e^"^ KHT)e^'^ K 


The input to the system is y^{t-T), the displayed variable to 
the operator, and the output is u(t), the controlled output of 
the operator. 

This model is based on the assumption that the operator has 
an "internal” aircraft model that is suitably modified by a 


PREC.vtJiNG PA 
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Kalman filter derived from the internal model and state noise, - 
An optimal controller is also implemented which ultimately 
produces the stick motion (output) of the operator. 

In the above equations, the following definitions are used 

(1) K = .Kalman gains = PH^ ' 

where P is the solution of the filter Riccati equation ‘ 

FP + pf"^ + - ph"^ hp = o 

T 

(2) X = optimal control gains = g"^ G S 

where S is the solution of the control Riccati equation 
SF + f'^S + Q - SG g'^ g’^S = 0 

(3) T = operator's "reaction” time (in seconds) 

The input to the algorithm are the matrices, covariances, 
etc. associated with the plant (e.g., F, G, H, W, etc.). At 
each iteration of the parameters to be estimated, F*, G* and 
X are computed as well as the sensitivity matrices 3F*/96, 
9G*/90, 3X/30 where 0 is a generic representation of the 
parameter(s) to be estimated. 

The sensitivity equation for F* is given by 

F* = (I + K'HT)(Fg-K'Hg- K'H - GXq - G^X) 

+ [TK'Hg + TK^H + TgK'H][F - K'H - GX] 

FT 

Where K' = e K = (I + FT)K = K + TFK 
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The sensitivity equation for G* is given by 


Also 


G* = (I + K'HT)(KJ) + (TK’Hg + TKJH + TgK’H)K’ 


= Kq + TFKg + TFqK + TgFK 


= (I + TF)Kg + CTFq + TgF)K 

W = s“ * s" «9S 

- g‘‘[G’'Sj H- gJs - gj|X] 


g-‘ gg g-‘ g’’s 


where Sg is obtained implicitly from 

(F - GA)’'Sg t Sg(F - GX) ’ - % ■ x’’^ggX - (Fg - GgX)’’S 


- sCFe - 


Depending on the parameter set, 6, to be estimated, the 
previous sensitivity equation(s) are constructed and computed 
in a sequence of various subroutines. 


In the SCIDNT implementation, F*, G* and A are computed 
transparently to the user. That is, at any iteration step, the 
matrices K', H, T, A are computed (as well as the sensitivity 
equation) and are used to compute F*, G* and A. Currently, 
the state equation matrices (and partials) are computed in 
subroutine STATG and measurement equation matrices (and partials) 
in subroutine MEASG. (The "general equation" mode of previous 
versions of SCIDNT are used for the generality required for this 
model.) 


STATG assumes that the user is using the general linear 
model (specified of the TITLE card). The four integers in 
columns 69-80 of the parameter cards indicate the position of 
each parameter in the system matrices. The parameters enter 
into the STATG routine in the P-array and the indices in the 

P-array to a set of local 
These matrices local to 

system matrix 
control matrix 
noise matrix 

process noise covariance 

measurement matrix 
measurement/control matrix 

measurement noise covariance 

weighting matrix (=Q) 

B (NQAC X NQAC) - control weighting matrix (=g) 

TAU (scalar) - operator's prediction/delay time 

(=T) 

IPDX elements are calculated by the user assuming the 
matrices are stacked columnwise in the order shown above. NQAC, 
Np!aC and NDAC are the number of controls, measurements and 
disturbances, respectively, for the plant model. F* , G*, X and 
DF, DG and DX (partials of F, G and X) with respect to the 
parameters) for the operator model are calculated by STATG from 
FAC, GAG, GAMAC, QAC, HAG, DAG, RAC, A, B and TAU. 

The entry points of STATG and their functions are described 
below: 


IPDX array act as pointers from the 
plant and human operator matrices. 
STATG are: 


FAC 

(NS X 

NS) 

- plant 

GAG 

(NS X 

NQAC) 

- plant 

GAMAC (NS 

X NDAC) 

- plant 

QAC 

(NDAC 

X NDAC) 

- plant 
(=W) 

MAC 

(NPAC 

X NS) 

- plant 

DAC 

(NPAC 

X NQAC) 

- plant 
(=0) 

RAC 

(NPAC 

X NPAC) 

- plant 
C = Vy) 

A (NS X NS) 

- state 
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STATG (main entry) - Calculates operator model F* and G* 

given FAC, GAC , GAl-lAC , QAC, HAC , 
DAC^(=0), RAC, A, B and TAU 

MEASG - Calculates X and D (=0) 

T 

GQGG - Calculates TQT for Operator model 

(since process noise is not imple- 
mented, this is always zero) 

it 

DERVG - Calculates F , G (operator) 

DGQGG - Calculates 8/86 (FQF^) for operator 

(always zero, since process noise is 
not implemented) 

DERM - Calculates X9 (operator) 

DERR - Calculates 8/80 for operator. 

Since IRCMP=0 is assumed, elements of 
Vjji are not identified. It is, 
however, computed "after the fact" by 
computing the standard deviation of 
the residual between the input control 
(stick output of the operator) read 
from data and the estimated control of 
the operator. 

Although restrictions on f TQT (Q is used 
internal to STATG, however this corresponds to external to 
STATG) have been pointed out above, the user can still identify 
elements in GAMAC, RAC, QAC existing in the previous versions of 
SCIDNT. 

All STATG entries interface with the existing SCIDNT 's 
UPDATE subroutine. Other SCIDNT’s subroutines changed from the 
previous version are: 

OUTERR - arrays depending upon number of controls, enlarged 
to accommodate 5 controls 

INREAD - reads controls/measurements 

In order to avoid overwriting intermediate results which 
might be in arrays DUM, DUH2, DDM in common block /DDM/ , all 
scratch storage used by STATG is in local array D. Dimensioning 
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information for D appears in the comments of the code for 
STATG. Other arrays which must be re-dimensioned (if number of 
measurements, states, controls, or process noise is increased 
above current maximum) are: 

FAC, DFaC (NS**2) 

GAC, DGAC (NS*NQAC). 

GAMAC, DGAMAC (NS*NN) 

QAC, DQAC (NN*NN) 

HAG, DHAC (NPAC^NS) 

DAG, DDAC (NPAC^NQAC) 

RAC, DRAG CNPAC*NPAC) 

AAG, DAAC (NS* NS) 

BAG, DBAC (NQAC*NQAC) 

P (NS*NS), K (MS*NPAG}, KPR (NS*NPAC) , 

IKPRHT (NS*NS), FKHGC (NS*NS) , 

DP, DK, DKPR (same as P, K, KPR), 

EVLRKF, EVLIKF, EVLROC, EVLIOC (NS) 

EVCKF, EVCOG, EVCIKF, EVCIOC (NS*NS) 

S, DS (NS*NS) C, DC (NQAC*NS) 

Also, matrices such as DUM, DUM2, DUMl , FAA, HAA, etc. must 
be re -dimensioned in 0UTER2RR. 

5.2 PROGRAM STRUCTURE 

In order to facilitate the future (potential) upgrading of 
the identification algoritlim, the following section describes the 
program structure generic to the original version of SCIDNT. 
Currently, the identification algorithm utilizes the output error 
method embedded in the subroutine OUTERR. Utilization of the 
subroutine DRIVER can be used to incorporate process noise in the 
model. 

The basic program structure is shown in Figure 5.1. 
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DRIVER and QUTERR are versions of the main routine. 
DRIVER should be used if process noise is present. DRIVER may 
also be used if process noise is not present, but in this case 
OUTERR may be used. OUTERR requires less storage than DRIVER 
to work with the same number of data points. 

Both routines read cards which specify page heading, and 
number of states, controls, measurements, and process noise 
sources. This information is passed to the subroutines so that 
arrays may be dynamically dimensioned. Currently, the maximum 
allowable values are 

number of states £ 5 

number of controls £ 3 

number of measurements £ 5 

number of process noise sources £ 1 

number of data points £ 301 for DRIVER 

£ 501 for OUTERR 

number of identifiable parameters £ 18 for OUTERR 

£ 10 for DRIVER 

S^'iAIN is called by DRIVER or OUTERR. SMAIN sets up the remaining 
program for identifying parameters. It reads cards which specify 
the parameter values, which parameters are to be fixed, bounds 
on the parameters, and other input data. SMAIN conducts the 
iterative maximum likelihood parameter search with successive 
calls to UPDATE. SMAIN also handles printout, printer plots, 
and writes time histories on tape or disc files. 

UPDATE is called by SMAIN. UPDATE solves linear ordinary 
differential equations for state time histories and for the 
sensitivity of the likelihood function with respect to parameters 
that are to be identified. It also computes the likelihood 
function itself. It calls STATG to determine state dynamics and 
control distribution matrices and their gradients in terms of 
the current estimates of the parameters. Similarly, it calls 
MEASG for the measurement equations. 
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STATG and MEASG are subroutines which define the matrices 
in the state and measurement equation, respectively. 

5.3 INPUT SPECIFICATIONS 

SCIDNT requires t\\ro classes of inputs. The first type, 
which will be referred to as "card- input , ” defines the number 
of states, controls, etc., and denotes the parameters of the 
model which are to be identified. SCIDNT reads card-input from 
unit 5. 

The second type of input consists of tabular values of the 
measurement and control time history. SCIDNT calls a subroutine 
INREAD once before beginning the identification algorithm. Sub- 
routine INREAD reads the values of measurements y and controls 
u for the entire time period of the experiment. The user must 
supply his own version of INREAD for "real" data. This is 
because in general, INREAD must read data in many different 
formats for various types of simulation or flight test data. 

For simulated data, INREAD is compatible with data written to 
disc from the TVOPT program. 

Card-input to the SCIDNT computer program consists of a 
sequence of cards constructed from various "card types." That 
is, several cards of a single type may be required in the input 
sequence. The card-input sequence is described in the following 
paragraphs. The card types are described in Table 5.1. 

Cl) One card of type 1. It contains run identification 
information which is printed at the top of each page 
of output. 

(2) One card of type 2. The number of states, measure- 
ments, controls, and process noise sources specified 
on‘ this card will be used in array dimension limits 
throughout the program. If they exceed the maximum 
dimensions allowed in the program, then a message will 
be printed and execution will stop. (The limits may 
be increased by re-coding DRIVER, the main routine.) 
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Table 5.1 

Input Card Formats 


CARD 

TYPE 

COLUMNS 

' 

FORMAT 

VARIABLE NAME 

DESCRIPTION 


1-76 

1 9A4 

TITLE 

Program identification informa- 
tion. Specify GENE as the first 
four characters. 

1 2 , 

1-5 

15 

KIMAX 

Maximum number of iterations 
(default value is 6) 

. ' ' 

6-10 

15 

K2MAX 

Maximum number of step cuts 
(default value is 4) 


11-20 

DIO.O 

DLl 

Sample time interval of time 
histories 

1 ! 

i i 

21-25 

15 

NNS 

Number of states (default value 
is 5) 

‘ 1 

26-30 

15 

NNP 

Number of measurements (default 
value is 5) 


31-35 

15 

NNQ 

I 

Number of controls (default 
value is 3 for a lateral case 
and 2 for a longitudinal case) 


36-40 

15 

NN6 

Number of process noise sources, 
must be >1 (default value is 1) 


41-45 

15 

NMAXP 

Total number of non-zero 
parameters 

3 

1-5 

15 

IPL0T 

, 

=3, if both printer plots and a 
tape of the variables T, Y, and 
YPL0T (time, measurements, and 
measurement estimates) are to be 
made 

=2, if only the tape is to be 
made 

=1, if only the printer plots 
are to be made 

=0, if neither is to be done 







m 
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Table 5.1 (Continued) 


CARD 

TYPE 

COLUMNS 

FORMAT 

VARIABLE NAME 

DESCRIPTION 


6-10 

15 

N0ISE 

=1, if process noise is to be 
included in the model 





=0, if no process noise is to be 
included in the model 


11-15 

15 

IRCMP 

=0, if estimate of measurement 
noise covariance (R matrix) is 
computed internally in program 
(usual procedure) 





=1, if estimate of measurement 
noise covariance is fixed to 
input values 


16-20 

15 1 

IPRNT 

=0, if extra printout for 
diagnostic purposes is not 
required 





=1, if extra printout for 
diagnostic purposes is desired 


21-25 

15 

IINF0 

=0, if the a priori information 
matrix is not to be input, but 
zero filled 





=1, if the a priori information 
matrix is input (see card type 7) 

1 j 

26-30 

15 

lEIGF 

=0, if the eigenvalues of F are 
not desired 





=1, if the eigenvalues of the 
initial and final F are to be 
computed and printed 


' 

31-35 

15 

IGUST 

=1, if estimate of wind gust is 
to be plotted (ignored if 
N0ISE =0) 






=0, if not 


36-40 


INC 

Printer plot interval 



1 
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COLUMNS 

FORMAT 

VARIABLE NAME 

DESCRIPTION 

1-5 

15 

IC0LK(i) 
i=l, NP 

=0, if the i^*^ column of K, the 
Kalman gain matrix, is to be 
zero filled 

20-25 

15 


=1, if the i^*^ column of K is to 
be retained (usual procedure) 

2-5 

14 

J1 

Parameter number as used in the 
model 

7-10 

A4 

i 

J2 : 

=blank, if the parameter is not 
to be identified 




=non-blank (such as an asterisk) 
if it is to be identified 

11-30 

020.4 

p(ji) 

Initial parameter value 

31-50 

020. -4 

PL(Jl) 

Lower bound on parameter 

51-68 

. 

018. 4 

PU(Jl) 

Upper bound on parameter. If 
both the lower and upper bounds 
are zero on input, then they 
are defaulted to -1000 and 
+1000, respectively 

69-80 

413 

IPX 

Parameter location 

1 

A1 

ECHK 

Any character (not a blank) to 
denote the end of cards of type 5 

1-80 

i 

801 0.10 

INFO(iJ), 
j=l ,. . .no. of 
parameters 
being identi- 
fied 

i row of the information matrix 
(continue the row on additional 
cards as needed). Repeat for 
each row of INFO, beginning each 
row on a new card. 

blank 




blank 





blank 




















If an input variable on this card has a standard 
value, then the variable will default to the standard 
value if the input yalue is zero (see Table 5.1). 

The sample time interval is data dependent and has no 
default value. The program computes its own integra- 
tion interval to strike a balance between computa- 
tional speed and accuracy (see section 5.4). 

(3) One card of type 3. This card flags various program 
options. 

(4) One card of type 4. This card specifies which columns 
of the Kalman gain matrix, K, are to be retained or 
filled with zeroes. It is infrequently but sometimes 
advantageous to delete columns of K to enhance 
convergence of the parameter search. The usual 
practice is to retain all the columsn of K. 

(5) Up to f#!AXP cards of type 5, one for each non- zero or 
identified parameter in the model. 

(6) One card of type 6 to designate the end of cards of 
type 5 in the input sequence. 

(7) If IINF0 = 1 on card type 3, then the a priori infor- 
mation matrix is input here by including as many cards 
of type 7 as are necessary. If IIMF0 = O', then no 
cards of type 7 should be used. 

(8) Three blank cards. 

• Card type 5 (see Table 5.1) is required for input of the 
system matrices. The matrices have the dimensions as shown in 
Table 5.2, 


Table 5.2 

Input System Matrices 


MATRIX 

DIMENSIONS 

F 

NS X NS 

G 

NS X NQ 

r 

NW X NG 

w 

NG X NG 

H 

HDx NS 

D 

NP X NQ 

h 

; NP X NP 

Q 

NS X NS 

g 

NQx NQ 

T 

scalar 


where NS = “ of states 

NP = # of measurements 

NG = # of process noise sources 

NQ = # of control inputs 

The user must allow for at least one column of T and one 
element of V^. 

The user defines the characteristics of both constant 
parameters and parameters to be identified using "card-input" of 
type 5. The last four integer entries define the location of 
the particular parameter in the system matrices. The location 
is an integer number specifying the vector location in the 
matrix list; i.e., each matrix is stored by column and the 
matrices are stored in the order shown in Table 5.2. Thus, 
location 1 implies F(l,l) and location NS*NS+1 implies G(l,l), 
etc. By using more than one of the four integer location inputs 
possible with each parameter, the user can direct SCIDNT to 
place a given parameter in up to four matrix elements. If a 



particular matrix location does not have a parameter value 
assigned, SCIDNT assumes a value of 0.0. 

In summary, the input required for each parameter according 
to the format defined above is: blank, parameter sequence 

number (1 to 38) , a non-blank character if this parameter is to 
be identified, the parameter value, allov;ed lower limit on 
parameter value, allowed upper limit on parameter value, and one 
to four location specification indices. 

An example parameter card format is given below: 

If 


NS = 5 
NP =» 5 
NQ = 2 
NG = 1 

then the parameter card 

1 * -2.62 -1000.0 1000.0 1 21 42 62 

will give these results for F and H: 


125 




[ 21 ] 


-2.62 0 

0 0 

0 0 

0 0 

0 0 


0 0 -2.62 

0 " 0 0 

0 0 0 

0 0 0 

0 0 0 



[42] 


[62] 


-2.62 

0 

0 

0 

0 


000 -2.62 
0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 


F(l,l) = location 1=1 
F(l,5)"’= location 5(5-1) + 1 = 21 

H (1 , 1) = location NS^ + NS • NQ + NS • NGN + NGN^ + 5(1-1) 
H(l,5) = location of H(l,l) + 5(5-1) =62 


2 = 42 


5,3.1 Subroutine INREAD 


The purpose of subroutine INREAD is to retrieve measurement 
and control time history data from cards, magnetic tape, or 
other mass storage devices. 

Calling Sequence: CALL INREAD (U, NQ, NDPMAX, Y, NP , T, 

TI, TF, NDPA) 

Inputs: NQ - number of controls 

NP “number of measurements 

NDPMAX * maximum number of data storage locations 
(i.e. data points) for each measurement 
or control 


Outputs : T 

= array for 
point 

U 

= array for 
histories 


storing the time of each data' 
storing the control time 


Y = array for storing measurement time 

histories 

NDPA = actual number of data points ready by 
the program 

Array Dimensions: T (NDPMAX) 

U(NU, NDPMAX) 


Y(NY, NDPMAX) 


Notes : 


(1) The program uses the array T only for the plot outputs 
The identification algorithm does not require T. 

(2) TI and TF are the initial and final times of the 
measurement time histories. The user can ignore these 
two arguments while writing INREAD, except that two 
dummy arguments must be included in the argument list. 

(3) INREAD should either read the value of NDPA before 
reading T, Y, and U or should count the number of 

• columns of these arrays as it reads them. The main 
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program does not initialize NDPA before calling 
INREAD. 



(4) 

INREAD should not reset the values of NQ, NP 
NDPMAX. 

, or 


(5) 

INREAD may perform any desired preprocessing 

U. 

of Y and 

5.4 

OUTPUT GUIDE AND EFFECTIVE USE OF THE PROGRAM 



Three forms of output can be obtained from Linear 

SCIDNT. 

They 

are : 




(1) 

Tabular printout showing the progress of the 
estimation process. 

parameter 


(2) 

Printer plots of time histories for measurements and 
measurement estimates using the final values of 
identified parameter and inputs. 


(3) 

Measurement, measurement estimates and input 
history data stored on magnetic tape for use 
generating off-line plots. 

time 

in 


The first type of output is always produced. Types 2 and 3 are 
optional and may be selected by the user (see input data cards 
description) . 


The option is provided to put measurement and control time 
history data on magnetic tape, primarily for the purpose of 
producing off-line plots analogous to the printer plots. Data 
is written to tape (currently designated as being on logical 
unit 2) in unformatted records as follows: 
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■;i 



RECORD 

DATA 



1 

NS, NQ, 

NP , NEW 



NS 

= Number of 

states 


NQ 

= Number of 

controls 


NP 

= Number of 

measurements 

; 

NEW 

= Number of 
variable 

data points per 


2 thru NEW+1 T(K), U(K), Y(K), YPL0T(K) 

T = Time for data step 
(1 variable per record) 

U = Controls for data step 
CNQ variables per record) 

Y = Actual measurements for k^^ 
data step (NP variables per 
record) 

YPL0T = Estimated measurements for 
k^“ data step (NP variables 
per record) 

NE¥+2 End-of-file mark 


The computer program's applicability is quite general. The 
following guidelines may make the program more useful for the 
evaluation of stability and control derivatives from test data. 
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5.4.1 Data Integration Step Size 

The program uses a second-order Runge-Kutta method to solve 
ordinary differential equations for the states and for the 
parameter sensitivites . This method is accurate only if the 
step size is much smaller than the smallest time constant of the 
system. However, a step size that is too small will require 
excessive execution time and may lead to problems with roundoff 
errors. A sampling rate of 20-50 times the highest system 
natural frequency (Hz) should be an adequate compromise. 

Linear SCIDNT automatically chooses the integration step 
size. Subroutine INTSTP sets the integration step size DTINT to 

DTINT = DTSMPL/K 

where DTSMPL is the sample period for the data and K is the 
smallest integer greater than ten times the magnitude of the 
largest eigenvalue of the system dynamics matrix F. This value 
is adjusted before each parameter identification iteration to 
reflect changes in the F matrix. 

I The user can make estimates of relative program running time 
using estimates of the largest norm of the system eigenvalues. 

5.4.2 Primary and Secondary Parameters 

The computer program can identify any or all unknown 
parameters if enough information is available about these 
parameters from the data. To identify a large number of param- 
eters, a very careful procedure is often required to ensure 
convergence. In the procedure which has been found to be most 
successful, the parameters are divided into two or more groups 
in the decreasing order of the effect they have on the response. 
Initially, only the most important parameters are identified. 
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leaving the remaining' ones fixed at a priori values. Once a 
reasonable convergence is achieved on these parameters, the 
second group may be added to the list of identifiable parameters. 

The identification is carried out using this new set of param- 

( 

eters, which are identified until a reasonable convergence is 
reached. This procedure is repeated until all the parameters 
are included in the identified set. 



5.4.3 Initial Parameter Estimates 

Good a priori values (e.g., from other \\?ind tunnel data or 
flight tests) should be used for start-up. If they are not 
available, it may be necessary to use a least-squares type of 
procedure to obtain starting values. 


5.4.4 Data Record Length 


Sufficient data length should be used to identify parameters. 
If the data length is too short, the identified model may yield 
a good time history match even though the parameter estimates are 
inaccurate. A data length equal to 2 to 3 times the longest 
period of the system should prove adequate. Keep in mind the 
required maximum number of data points. 




5.4.5 Data Sampling Interval 


In order to assure adquate information content in the data, 
a sampling interval of at least 25 times the highest system 
natural frequency is required. A faster sampling rate provides 
somewhat more accurate parameter estimates because of some addi- 
tional information available in the sampled data. However, the 
algorithm realizes diminishing returns in terms of increased 
parameter accuracy as the parameter estimates approach the 
continuous time case. 
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5.4.6 Diagnostics 


ABNORMAL TERMINATION 

CORRECTIVE ACTION 

CONDITION 


i 

Step Cut Limit Exceeded 

CD 

Increase step cut limit and 
restart using parameter estimates 
from final iteration of terminated 



run. 

. 

(2) 

Delete parameters with low confi* 
dence (F-value small) from the 
search. 


Iteration Limit Exceeded (1) Restart program using parameter 

values from final iteration of 
terminated run. 

C2) Increase iteration limit. 

NORMAL TERMATION, BUT CORRECTIVE ACTION 

PROBLEMS IN RESULT 


High Value of a Partic- (1) 
ular Measurement's 
Noise Covariance 


Faulty data: Inspect for correct 

sign, slipping, lost values, 
"shot” noise. 


High Standard Deviation 
on a Particular 
Parameter Estimate 


(2) Incorrect values for "fixed" 
parameter . 

(1) Insufficient excitation of the 
mode ivhich that parameter 
influences. Rerun identifying 
only the more important parameters 


(2) Too many parameters identified 
for that particular data record 
Ci.e., overparameterization). 
Rerun, identifying only the m.ore 
important parameters. 

C3) Parameter influence masked by a 
parallel influence from another, 
much larger parameter (i-e*» 
identifiability) . 
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5.5 PROGRAM FLOWCHART (STATG) 


The program flowchart for STATG is shown in Figure 5.2. 
Note that STATG uses the VASP routines 


TRANP 

UNITY 

INV 


These must be present in a user library when the program is 
loaded and executed. Also, STATG usies the extended VASP routines 


ALA 

APALB 

APALBT 

ATXB 

AXB 

AXBT 

AXNXBT 

AZERO 


j These routines were used because they have more general 
calling sequences, and, in many cases, combine the operation of 
transposition with multiplication or addition, resulting in more 
efficient use of temporary storage and fewer subroutine calls 
for complicated matrix algebra expressions. Descriptions and 
calling sequences for these routines are as follows (N followed 
by a matrix name, e.g. NA, denotes a 2-element integer array 
containing the numbers of rows and columns of the matrix) . 


CALL ALA (A, NA, B, ALPHA) (B = ALPHA*A) 

Scales A by ALPHA, storing result in B. A and B can be the 
same matrix. 
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"a" subscripts 
denote FAC, etc. 
T = TAU 



Figure 5.2 Program Flowchart CSTATG) 


F,P + 

a 


CALCULATE P AND K MATRICES 
(INNER) 

« T T_-i 

r-.Q.r‘ 


pfI + 

a 


a^a a 




K = 


T -1 

ph‘r 

a 




= 0 
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Figure 5.2. Program Flowchart (STATG) (Continued) 











SL. 



Figure 5.2 Program Flowchart (SJATG) CConcluded) 










CALL APALB (A, NA, B, C, ALPHA) (C = A + ALPHA*B) 

Scales B by ALPHA and adds it to A, storing result in C. 

A, B, and/or C may be same matrix. Used for matrix addi- 
tion and subtraction when ALPHA = 1 or -1. (No multiplies 
are performed in these cases.) 

CALL APALBT (A, NA, B, C, ALPHA) (C = A + ALPHA*B^) 

Scales transpose of B and adds to A, storing result in C. 

A, B, and/or C may be same matrix. Note that A + B must 
be square . 

CALL ATXB (A, NA, B, NB, C, NC) (C = a”^*B) 

Multiplies transpose of A by B, storing result in C. A + C 
or B + C may not be same matrix. NC is output parameter. 

CALL AXB (A, NA, B, NB, C, NC) (C = A*B) 

Multiplies A by B, storing result in C. A + C or B + C may 
not be same matrix. NC is output parameter. 

•; ! 

CALL AXBT (A, NA, B, NB, C, NC) (C = A*B^) 

Multiplies A by transpose of B, storing result in C. NC is . 
output parameter. A + C or B + C may not be the same 
matrix. 

CALL AXWXBT (A, NA, W, B, NB, C, D) (C = A*W*B^) 

Multiplies A by W by transpose of B, storing result in C. 

D is an NA(2) x 1 scratch vector. W is assumed to be 
dimensioned (NA(2) x NB(2). A + C Or W + C or B + C must 
not be same matrix. D must also be a separate storage area 
from A, W, B, or C. 

CALL AZERO LA, NA) (A =0). 

Zeros out matrix A. 


f 


I 

f 

.1 

( 


I 


n 

i '' 
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5.6 TSS/360 OPERATION 

I . ' 

The following control cards are required to execute SCIDNT 
on the NASA ARC TSS/360 system. 

(1) L0G0N... 

(2) DDEF FTlOFOOl, VS, TIME .HISTORY, DISP=0LD 

(3) JBLB SCIDNT 

(4) JBLB HUMAN0P 

(5) L0AD BL0CKD$$ 

(6) CALL 0UTER2$$ 

(7) {input data 
('8) L0G0FF 

The description of the above execution deck follows: 

(1) : Usual L0G0N card. 

(2) Data definition card specifying the file on which time 
history is stored. 

(3) ! Assign library SCIDNT. 

(4) Assign library HUMAN0P. 

(5) Load in block data. 

(6) Execute program. 

(7) Data. 

(8) Usual L0G0FF. 
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APPENDIX A 


The following program is a "stand-alone” computer program 
which can be used to create a general output disc file to be read 
as input for the TVSR program. Specifically, the program reads 
the input data from cards and writes to disc file in the 
identical format as does the TF program. Since the TVSR program 
inputs data in identically this format, this stand-alone code can 
be used to generate input data for TVSR without the need to 
exercise the TF program. 

' The input format for TVSR is given in Table 3.2. The 
folloi\ring program can serve to create this data from an input 
card deck; an example is also shown. Note: In the example, 

card 1 indicates 3 x's. These x's correspond to (1) the number 
of plant states, (3) the number of human states and (3) the 
niimber of augmented states. In the case when the input vector 
is not augmented (with the deterministic input, p(t)), these x 

i 

values are arbitrary. For the case when the input vector is 
augmented, care should be taken in inserting the correct values. 


r»oo r> ri n nnn nr»o n n r* 


ImPLICM Pf-AL*H IA-h»0-2i 
DIMFNSION IFKT (3) 

DlMKNSlOF! Ftl^H5)i r.MSt|5), VyfStS) 

dimension VA(‘i|SJ 

DATA IfhT / '(SEt, ilO.S'f •)» / 

READ OIMEMSIOkS 

REAOrStlOOl NS< |INP» NSHf NMHi NSCt NJt NHt NSA« NCA 
too FORMAlOI'il 

head matrices at siahi time 

call HDMAT(F,NS,NS, IFMT,5) 

Call homa1|G,ns,ni , ifmt.S) 

Call rdmat(H.nm,ns, ifmt.si 
Call hom a I ( v. , nnp # nnp , i fht , 5 j 
Call RDMAlcVy.NMHtNMH»IFMTi5) 

Call HOMAlf VA» NCAt NCA# IFHT# 5 ) 

WRITE UjMENSInNS AND HatRICFS AT START TIME 

WRllE(lO) kS, MNP, msH# NhH, NSc# NI# NH, nSa# NCa 
Call BKTMATiFfNSfNS, |0) 
call BWIMAt(r,#NS#NI,lO) 

Call HWTMAlCh.NM.NS, iOT 
call pWlMAl(W#NNP»NMP» 10) 

Call HWIMAKVYtKMU.NMH# IO) 
call RWIMAKVA# NCA, NCAt 10) 

Endfile Jo 

read matrices at end Time 

Call RDMATiF.NS.NS. IFt1T*5l 
Call RDMAlfG.NStMI# I fHT»5) 

Call ROMAI(h,NM,NS,IFMT#S) 
call RDMA1(W,NNP|NNP,IFMT,5) 

Call RDMAI(V¥,NMHfNMH»IFMT,5) 

Call rDmau va. nca, nca, ifmt# "j) 

Write dimensions and matrices at end HMe 

WRIlEdlJ N'S, nnp, NSH, NmH, N'Sc, NI, NM, N3a, NCA 
call hN|MAi(F.N3,NS,| 1 ) 

Cam 8l'TMAT(r.,NS,NI,M ) 

Call bwimai(h,nm»ns, i i j 
CALL nHTHAt (W,NNP»NNP# 11 ) 

Call 0WlMATfVT#NMH,NMH, M) 

Call rwtmatIva# nca, nca, iti 

FNOFILE It 

STOP 

End 
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KOU BC - ( * Hj; 

iMPLICit BEAl-^e 

Dimension *(nk,nci 

Do too UIiNR 

READCIUNIT.IFMI) (A(l*J)tJ=t,NC) 

too Continue 
Return 
End 



CZIII? | ‘” ^ '* " * * ** '| gww^ 



O 


I— 

I - j c, -_ 

o 

03 
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SUBROuUNt HWTf<Ar(AtNH»Nc*lUNlT) 
iMPUcn REAL*a(A-M|0«2) 

DIMENSION a(NRiKC) 

MRlTEdUNJt ) r(A(I.JJtU|rNfi)f J = l>NCj 
I Return 

! END 

I 






I, — - 




11 


X i 1 X 

1 7 

5 5 

X 3 i 



2 0.0 

0.0 

0.0 

0.0 

0.0 

0.990 

0.0 

3 0.0 

0.0 

- 0.263 

0.0 

0.0 

0.0 

1.0 

4 0.0 

0.263 

0.0 

0.0 

0.0 

0.139 

0.0 

5 0.0 

0.0365 

0.0 

0.0 

0.0 

1.0 

0.0 

6 0.0 

0.0 

- 0.263 

0.0 

0.0 

0.0 

1.0 

7 








8 








9 

10 


7 . E - 70.0 

0.0 

0.0 

0.139 




7 . 

E -7 





11 p.o 

0.0 

0.0 

0.0 

- 0.990 




n 2 . 88 E -4 

t 3 • - 2 . 88 E -4 

14 , 

15 i 

16 1.0 - 1.0 

17 ^ 1.0 - 0.1405 1.0098 

18 1.0098 - 0.1405 

19 -.707 -.707 

20 .707 -.707 

21 1.0 
221.0 

23 

24 

25 
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